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SECTION I 


INTRODUCTION 


The mathematical operation of convolution is a very important function in 
many of the pattern recognition or pattern classification algorithms that are 
used to synthesize useful information about the earth from the electrical out- 
puts of satellite-borne environmental sensors. Since charge coupled devices 
(CCDs) are so well suited to performing convolution, it is not surprising that 
substantial Interest has arisen within NASA in developing CCDs for preprocess- 
ing environmental sensor data preparatory to Its transmittal to the ground. 


The one complication that limits the immediate usefulness of CCDs is the 
fact that the satellite preprocessors are learning machines and must adapt to 
the terrain they are observing. In terms of CCD hardware, this means that the 
convolutions must be programmable; i.e., the impulse responses of the filters 
must be electronically variable. 

The canonical CCD transversal filter, illustrated in Figure 1, has fixed 

weighting coefficients h^ n = 0, N - 1 that are coded into the photomask with 

which the device is made. The impulse response is determined by the h^ and is 

not electrically alterable. Various approaches have been tried to construct 

2 3 

filters with electronically alterable h^, ’ but all such attempts have failed 
to maintain the electronic simplicity of the fixed weighting coefficient trans- 
verse 1 filter. 


This program addresses two different techniques for performing programmable 
correlation 

• Frequency domain correlation using the CCD chirp z-transform (CZT) . 

• Time domain correlation using direct analog/analog correlation and 
binary/ana 1 og correlation. 

Two separate ICs have been developed under the program: 


32-point CCD CZT 1C. 



iagram of a Transve 
eighting Coefficien 
















32-point CCD analog/analog correlator. 


These two ICs are discussed in Sections II and III, respectively. In addition, 
a prototype IC capable of correlating a digital reference word with an analog 
input signal was evaluated under this program. The evaluation of the binary/ 
analog correlator is discussed in Section IV. 

A . Frequency Domain Correlation 

This approach to programmable correlation is based on the fact that corre- 
lation in the time domain is equivalent to multiplication In the frequency 
domain. Consequently, the frequency domain correlator Is based on taking the 
discrete Fourier transform (DFT) of the input signal. 

The discrete Fourier transform operation can be performed using the chirp 

4 

z-transform algorithm. The CZT gets its name from the fact that it can be 
implemented in an analog manner by (1) premultiplying the time signal with a 
chirp (linear FM) waveform, (2) filtering in a chirp convolution filter, and 
(3) pos tmu 1 t i ply i ng with a chirp waveform. 

When implemented digitally, the CZT has no advantages over the conventional 

4 

fast Fourier transform algorithm. However, the algorithm lends itself natu- 
rally to implementation with CTD transfersal fi Iters. 


Starting with the definition of the DFT 


F 


k 


^ -i2rrnk/N 
L f e 

n ^ 
n=0 


and using the substitution 


2nk = n^ + k^ - (n - k)^ , 


the following equation results 




• i2,., ,N-I 

- irrk /N 


S (f giTr(k-n)^/NJ 


“-n^O 


(0 


( 2 ) 


( 3 ) 
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This equation has been factored to emphasize the three operations that make 
up the CZT algorithm. It is illustrated in Figure 2. 


To implement the conventional N-point CZT, the CCD filters are chirp 
filters of length 2N-1 that chirp from -f^ to and the premultiply wave- 

form has a time duration N/f^ and chirps from zero to -f^. A physical inter- 
pretation in terms of correlation of the input chirp with the filter is given 
in Figure 3. When the input signal has zero frequency, the product with the 
premultiply chirp results in an input waveform to the filter that chirps from 
0 to -f^. The samples corresponding to frequencies near f = 0 are clocked 

into the filter first, and those near f = -f are clocked in last. This 

c 

sequence of samples results in a correlation peak at t = t^, when the product 
waveform has been clocked into the first half of the filter. When the input 
frequency is f^ ^ 0, the product with the premultiply chirp results in an 
input to the filter that chirps from f^ to -f^ + f ^ . The input waveform 
(V.^ X chirp) in Figure 3 corresponds to an input signal at a frequency f^ at 
time t = t^. This waveform is shifted to the right as t increases, resulting 
in a correlation peak at t^ The shift in time relative to the dc correlation 
peak is 



(^) 


In this way, the time axis of the output is calibrated in frequency. 


The CZT is very convenient for performing programmable correlation. 

Figure 4(a) shows that correlation of x(t) with h (t) can be accomplished by 

(1) performing the discrete Fourier transform (DFT) of x(t) to obtain X(f); 

(2) multiplying X(f) by H(f), where H(f) is the DFT of h (t) ; and (3) performing 
the inverse DFT on Y(f) = X(f) • H(f). 


4 
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Figure 4(b) is obtained from 4(a) by replacing the DFT and inverse DFT 
blocks with their CCD CZT equivalents, and 4(c) is obtained by simplifying 
4(b). Figure 4(c), therefore, represents an economical and useful way to 
perform programmable correlation in pattern classification machines, and this 
program had as one of its goals the development of CCD CZT ICs to perform 
spectral analysis for a broad range of NASA requirements including programmable 
cor re 1 a t i on . 

The 32-point CCD CZT 1C developed under this program is discussed in 
Section II. 

B . Time Domain Correlation 

An alternate approach to achieving programmable correlation is the direct 
time domain correlator. The analog/analog correlator is illustrated in Figure 5.^ 
It is similar in operation to the transversal filter illustrated in Figure 1, ex- 
cept that the weighting coefficients are electronically programmable by clocking 
an analog reference signal into the device. It operates by first loading CCDl 
with a reference signal and then clocking a second signal through CCD2. Since 
the charge is stored dynamically in CCDl, it decays due to thermal leakage and 
must be refreshed after 10 ms to 100 ms. 


Each stage of the CCD is tapped with a diffusion whose potential takes on 
a value approximately proportional to the signal charge, and four-quadrant 
multiplication is accomplished using the circuit of Figure 6. Transistors 
and act as voltage-controlled resistors, and the balanced configuration 
overcomes the inherent nonlinearity of MOS transistors used as resistors. 


The output V 2 from CCD2 Is applied to the gate of Q^, and the gate of is 

biased to V which is equal to the value of \l^ when CCD.- has the zero level 
ref 2 2 

of charge (half-full well). The output from CCDl is buffered by a source- 


follower whose output M ^ is applied to the drain of and Q^. The sources of 
and 0^ are connected to virtual ground current buses L and S'** whose 
currents are differenced in an off-chip differential amplifier. 


8 













V 


DD 



Figure 6. -Circuit Schematic of the Four-Quadrant Multiplier Used in the Analog/ 
Analog Correlator of Figure 5* 
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The currents flowing in and 0^^ are 


I3 =P(Vr,f - \-i\) \ > 


( 5 ) 


and 


li, = P(V2 - 


( 6 ) 


The difference current 


!. - I, = P(V„ - y ) V 

4 3 2 ref a 


(7) 


contains no nonlinear terms. 

The 64-point ana 1 og/ana log time domain correlator designed, fabricated, 
and tested under this contract is discussed in detail in Section III. 


An alternate approach to the analog/analog correlator is the binary/analog 
correlator. The concept of the binary/analog correlator is illustrated in Figure 
The weighting coefficients (h^h. . . , h^ ^ of Figure 1) are made electronically 
programmable by decomposing each one into a binary representation, h^ is repre- 
sented with M-bit accuracy by 


k -k 

h = T. h 2 ^ 

^ k=0 


( 8 ) 


For M“bit accuracy, M parallel CCD binary/analog correlators are put on a chip 

as shown in Figure 7- The most significant bit h^ of each coefficient is loaded 

into the static shift register (coefficient store) shown as elongated rectangles 

in the filter at the top of the figure. The second most significant bit h^ of 

each coefficient is loaded into the second coefficient store, and the least 
M-1 

significant bit h is loaded into the coefficient store shown at the bottom 
n 

of the figure. The analog input signal to be filtered is applied without 


Out 
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Figure 7.- Block Diagram of the Binary/Analog Correlator, 













attenuation to the top filter (most significant bit). The input is attenuated 

by a factor of two at the input of the second filter (second most significant 

H- 1 

bit), and it is attenuated by a factor of 2 at the input to the bottom filter 
(least significant bit). This attenuation is performed using capacitive ratio 
techniques similar to those employed in MOS MDAC technology. The coefficients 
stored in the static shift registers then control the weighting of the trans- 
versal filters, and when the outputs of each filter are summed together as 
shown, the result is 


H(z) 


V (z) 
out 

V. (z) 
in 


N=1 




N-1 

Z 

h° Z 

-n 
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n=0 
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(h'^ 

2~k^ 

)] ^ 

n=0 

-k=0 

\ n 



N-1 





Z 

h z 

-n 



n=0 
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+ 2 


-(M-I) ,M-I ^-n 

n^O 


( 10 ) 


( 11 ) 


The binary/analog correlator is discussed in detail in Section IV. 


13 



SECTION II 


32-POINT CCD CZT 1C 

This section describes in detail the design and operation of the frequency 
domain correlator. The design of each component of the CZT 1C is discussed, 
and applications and experimental results are described. 

A. Chip Overview 

The frequency domain correlator contains all the hardware needed to 
implement a 32-point CZT, as discussed in the introduction and shown schemati- 
cally in Figure 4. The 1C, shown in block diagram form in Figure 8, is com- 
posed of the following: 

• ROM 

• Up/down counter 

• Four MDACs 

• Four operational amplifiers 

• Eight output buffers 

• Clock generators 

• Four 64-stage CCD filters. 

A photomicrograph of the 1C is shown in Figure 9. 

2 2 

The 1C measures 6.04 x 5.69 rnm (238 x 224 mils ) and was fabricated 
using Tl's two-phase silicon gate coplanar electrode structure, shown in 
Figure 10. To provide storage and directionality, n-type ion implants (wells) 
are diffused under the two clock electrodes. As seen in Figure 10, when both 
the 0^ and 0^ clocks are off, charge (signal) resides in the wells (under the 
0^ electrode for this example). When 0^ is turned on, a potential energy 
gradient is formed, forcing the charge to move under the electrode. Next, 
both clocks are turned off, and the charge remains in the 0^ well. When 0^ is 
turned on, the charge moves under the 0^ electrode and has moved one bit down 
the CCD structure. 


14 




MDAC 






















17 


Because of the complexity of the analog circuits surrounding the CCD 
filters, much circuit modeling and MOSFET modeling was needed. To obtain 
accurate MOSFET models for use in the circuit modeling, data were taken on 
MOSFETs from previously processed ICs. These data were tabulated and fitted 
to the Sh i chman-Hodges MOS equations using the computer, and parameter values 
were extracted for circuit modeling. This procedure allowed the design of 
the peripheral circuitry to be accomplished with a high degree of confidence. 

In an effort to provide the greatest possible flexibility, the 1C was 
partitioned into three main signal processing sections. The first two sections 
are both dual MDACs. The MDACs and their associated circuitry are shown in 
Figures 11 and 12. The ROM is common to both sets of MDACs and has 16 output 
lines, 8 for the sin chirp words and 8 for the cosine chirp words. The sin 
chirp is available to one MDAC in each figure, and the cos chirp is available 
to the other two MDACs. The eight-bit external word is available to all four 
MDACs and may be individually selected for each set of dual MDACs. Each MDAC 
consists of an input sample-and-hold buffer, ROM sense and MDAC drive circuitry, 
the MDAC capacitor array, and an output buffer. The third section of the par- 
tition, shown in Figure 13, consists of two CCD input amplifiers (A^ and A^ ) , 

CCD input blanking circuitry, four 64-stage CCD transversal filters, two 
differential current integrators (CCD output amplifiers), and four output buffers. 
The CCD input amplifiers make the inverse of the signal available to the CCDs 
so that both CCD outputs to the DC I may be summed the same way. The CCD input 
blanking circuitry uses the T^ and T^ pulses from the up/down counter to load 
"zeros" into the CCD for appropriate periods. The T 2 pulse is always used. 

The T^ pulse is used only to implement the noncircular correlation function 
and must be stitch-bonded on-chip if it is required. 

Each of the three partitioned sections utilizes output buffering circuits 
so that the sections can be connected in several configurations using one or 
more ICs. Also, each section has the circuitry required to process real and 
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Figure lU’-Block Diagram of Internal ROM, Up/Down Counter, and Dual MDACs. 
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Figure 12. -Block Diagram of Second Dual MDAC . 
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Figure 13. “Block Diagram of CCD Input Amplifiers, CCDs, dCI. and Output Buffers. 






imaginary signal components. In addition to the three partitioned sections 
described above, the chip contains clock-generating circuitry. The clock 
generators derive all the internal timing pulses from the two-phase clocks 
input to the 1C. 


B. CCD Filter 

The four filters needed to implement the CZT algorithm are chirp con- 
volution filters with 2N-1 stages. A photograph of the four CCD filters is 
shown in Figure l4. The weighting coefficients for the sine and cosine chirp 
f i 1 ters are 


h 


cos 

k 


cos 



k = 0,31 


h 


s i n 
k 


s i n 



k = 0,31 . 


The coefficients were realized using the split-electrode technique, and 
because the CZT !C was designed to perform the true CZT, the filter tapweights 
were repeated for the remaining N-1 stages of the filters. 


The CCD filter structure is Illustrated in Figure 15. Each filter is 
surface channel and uses two-phase clocking. Each stage is 25.4 pm (1 mil) 
long and consists of a 12.7 pm (0*5 mil) long first poly electrode and a second 
poly electrode that is effectively 12.7 M-^ (0.5 mil) long. The barriers are 
5.1 pm (0.2 mil) each, and the wells are 7-6 pm (0.3 mil) each. The channel 
width Is 155 pm (6.1 mils). The split electrodes are the first poly electrodes. 
The split electrode gaps are 5 pm (0.2 mil), and the channel stop beneath each 
gap is 12.7 M-m (0.5 mil) wide. 

The channel stop beneath each gap serves two purposes. The first is to 
prevent the electrons of the signal charge packets from being trapped in these 
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Figure I^. • Photomicrograph of CCD Transversal Filte 















areas and causing serious degradation of the charge transfer efficiency. The 
second reason for the channel stop is that the position of the weighting 
coefficients is defined by the same photomask that defines the CCD channel. 
Thus, the weighting coefficient accuracy is determined by the accuracy with 
which the photomask is made and is not dependent on alignment of different 
photomask levels in the device fabrication process. With this approach 
approximately 0.16% accuracy in the weighting coefficients is possible with 
a 155 M-m (6,1 mil) wide channel in the CCD. This gives an area for the first 
well area of 7-6 x (1.55 ” 12.7) = 1084 pm^ [0.3 x (6.1 - O. 5 ) = 1.68 mil^]. 

The CCD clocking sequence is illustrated in Figure 16 relative to the two- 

phase input clocks and • The analog input to the CCD is capacitively 

coupled to the input diode. The fat-zero level on the diode is set to 

when 0| comes up. After 0| falls, the CCD input signal is allowed to change 

and the IP/G signal pulls up, dumping signal charge into the receiving well, 

2 

which is controlled by a dc bias. The total receiving well area is 7 IO pm 
2 

(1.1 mil ), making the capacity of this well 65 % as large as the capacity of 
the following wells. This prevents saturation of the filter. The signal is 
sampled when IP/G falls. Next, the CCD first 0^ clock rises and falls, dump- 
ing the sampled signal charge under the first split electrode. The CCD first 

+ 

0^ electrode rises to a level less than that of the CCD 0^ clock lines to 
make sure that all charge is transferred. At the same time, the CCD 0^ 
electrodes fall, transferring the charge packets under the CCD 0^ electrodes 
to the remaining split electrodes. These split electrodes must be floating 
when this occurs in order to sense the charge packets. 

+ 

The capacitive coupling of the CCD 0^ line to the CCD 0^ lines is com- 
pensated for by another clock, CCD0^ , which is capacitively coupled to the 
CCD 02 lines by approximately the same amount. The CCD0^ clock comes up when 
CCD 0^ clock falls. The output signal is sampled before the CCD 0^ line rises 
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again. When the CCD 0~ electrodes are pulled down, the charge packets are trans 
ferred to the following CCD 0j electrodes. When each charge packet reaches the 
end of the CCD, it is dumped onto the output diode, which is tied to An 

effort was made in laying out the four CCDs and their interconnects to match the 
capacitance on all four output lines. Consideration was given both to capaci- 
tance to the substrate and to capacitance to the other clock lines. 

C . Ope ra t i ona 1 Am p 1 i f i e rs 

1 . Des i qn 

In the realization of the 32“point CZT, four operational amplifiers 
are needed to perform two functions. The first function is to provide gain 
and the inverse of the signal to the inputs of the CCDs. These amplifiers 
are called the input diffamps and are designated as Al and A2 in Figure I 3 . 

The need for the inverse of the signal can be understood best by considering 
the function of the differential current integrator (DC I ) . 

The general DC I scheme and the capacitances associated with the CCD 
are shown in Figure 17. The following capacitances are indicated in the 
f i gu re : 

+ 

C_ = stray to substrate, made equal on V nodes. 

^ ± 

Cqp = fixed overlaps from phase 1 (0^ to V nodes, also made equal. 

C (1 ± h.)/2 = overlap to 0^ for each stage, depends on h.. 

Cqp + N/2 Cqi^ = a coupling to 0^ equal for each node V . 

C (1 ± h.)/2 = \l~ electrode to channel capacitance for each stage. 

Ua I 

^D ~ = depletion layer capacitance for each stage. 

A possible timing scheme is shown in Figure 18. The operation is as 

+ 

follows. The split electrode nodes V are clocked to the voltage via the 
series switches while the phase 1 (0^) clock is high. The switches are 

then turned off, and nodes V~ are left floating. The next point in the timing 
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diagram is the turn-off of phase 1. At this point the charge transfers under 

+ 

the split electrodes. The nodes V are isolated, and feedthrough from (f), via 

+ 

the overlap capacitance could reduce the voltage V ^ \f 2 to a value low enough 
so that incomplete transfer could take place. To prevent this, a capacitive 
coupling to 0j is provided, which is approximately equal to the coupling to 0^ 

As the signal charge Q_ flows under the split electrodes, the node 
voltages V change in proportion to the sum of the tap weights (1 ± h.). The 
node voltages then are 


="2 - \l^0F " I/' "h') 

whe re 

C* = C 5 + 2 Cq|- + Cql j (2 + h) + Cg I {) ± h) , ( 13 ) 

and in series, and 

N 

h = S h./N, 
i=l ' 

= bias voltage for zero signal (fat zero), 

V. = signal voltage sampler, 

+ 

The node voltages V are to be sensed by the differential feedback amplifier A 
when the switches SA are turned on at the next stage in the timing diagram. 

The negative feedback of the amplifier causes its output to change until the 
input voltages and V are equal. 
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The output voltage 's proportional to the difference: 

- V- = Nh Cg^/C + Nh Cq/C + h,V,/C , (14) 

where C and C have been made equal by the addition of (C^. + C^) Nh to node 

+ 0L_ E" 

V . In practice, the addition may have to be made to V if h is negative. 

Only the third term In Equation (l4) is the desired output of the 
convolution filter. The first two terms need to be eliminated. This can be 
accomplished in two ways. The first involves an additional channel in parallel 
with the signal channel. This new channel is of width jh] • if is the 
signal channel width, and is located on the side of the V electrode if h >0. 
The new channel only carries a signal equal to V^. Again, in practice, two 
such channels may be required, one on each side of the main signal channel. 

Their widths are then adjusted so that they differ by jh] W^. A possible lay- 
out of such a structure is shown in Figure I 9 . The addition not only cancels 
the two unwanted parts in the output [Equation (14)], but also adds just the 
required capacitances to equalize C and C . 

A second way the two additional terms in Equation (l4) can be cancelled 
is appropriate to the chirp z-transform. Here the output of a cos and a sin 
chirp are added. However, the values of h for the sin and cos chirp are equal, 
at least for an even number of stages. Even when apod izat ion is added to the 
CCD chirp filter, as may be done for a sliding transform, the values of h 
differ by only a small amount. For example, for even N, h(cos) = h(sin) = 

1/^^. For a Hamming weighted filter we calculate 
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Figure 1 9.- Symmet r i zed and Compensated Two-Phase CCD Layout, 



N 

h (cos) 

h (s in) 

h (cos 

- sin) 

512 

0.0313 

0.0311 

1.7 X 

10“^ 

256 

0.0444 

0.0439 

5.0 X 

10“^ 

128 

0.0632 

0.0618 

1 .4 X 

io"^ 

64 

0.0903 

0.0863 

4.0 X 

10“^ 


0. 130 



-2 

32 

0. 119 

1 . 1 X 

10 

16 

0. 189 

0. 158 

3.1 X 

10-2 

8 

0.274 

0. 193 

8. 1 X 

10"2 


The difference is small for the larger transforms. 

The required equalization of capacitances and cancellation of 
unwanted signals can then be obtained by substracting the output of a cos and 
sin chirp filter. In practice, this subtraction is achieved by cross-connecting 
the v"*" (cos) and V (sin) nodes. No additional channel is needed. However, 
one of the input signals has to be available in both normal and inverted form. 
Since the length of our transform was small and gain was needed at the inputs 
to the CCD filters, we chose to use the second approach. Two operational 
amplifiers are used at the input to the filters and two amplifiers are used 
as DC Is at the output of the filters. The designs of the two diffamps are 
identical except for the values of the compensation capacitors and the 
feedback capacitors. 

The basic amplifier is shown in Figure 20 with the design parameters 
listed in Table I. A photomicrograph of the amplifier is shown in Figure 21. 

The amplifier has two stages of differential gain. Each stage is 
cascaded and uses depletion loads. Between the first and second stages are 
buffers (Ml4 and M15) which also serve as level translators. After the second 
gain stage, the outputs are doub le-buf fered by M30 and M31, and then by M33 
and M34. The current source biases of all buffers are cascaded to provide a 
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Figure 20 ." DCI Amplifier Compensated for Gain of 30 . 





















Figure 21.- Photomi crograph of Operation Amplifier Used on CZT 1C. 









Table I 


Vjjp = +20 at 3.5 mA 
'^SS " ° 

FET Geometries and Supplies 

Vbb=-5 at 1.5 

Cl = C2 = 0.075 

FETS 

W/L (mils) 

W/L (pm) 

Ml , M7, M23 

0,35/0.3 

8. 9/7-6 

M2 

1.25/0.3 

32/7.6 

M3, M4 

18/0.3 

457/7.6 

M5, M6 

1.30/0.3 

33/7.6 

m 8, M9 

0.35/0.7 

8.9/17.8 

M10,M1 1 

0.7/0. 3 

17.8/7.6 

M12, MI3 

0.75/0.5 

19/12.7 

M14, MI5 

1.5/0. 3 

38/7.6 

M16, MI7, M21, M22 

0.45/0.5 

11.4/12.7 

M18 

0.6/0. 3 

15.2/7.6 

MI9, M20 

2.45/0.3 

62/7.6 

M24, M25 

0.35/0.5 

8.9/12.7 

M26, M2 7 

0.5/0. 3 

12.7/7.6 

M28, M29 

0.55/0.5 

14/12.7 

M30, M31 

0.95/0.3 

24/7.6 

M32 

0.7/0. 5 

17.8/12.7 

M33 , m 41 , m 42 

0.65/0.5 

16.5/12.7 

M34, M37 

0.35/0.65 

8.9/16.5 

M35 

0.75/1.0 

19/25.4 

M36 

0.65/1.7 

16.5/43 

M38 

0. 8/0.4 

20/10 

M3 9, M40 

0.65/0.3 

16.5/7.6 

m 43 , m 44 

5. 1/0.3 

130/7.6 

m 45 

0.75/1.7 

19/43 

M53 

0.45/0.3 

11.4/7.6 

M54, M55 

0.35/0.85 

8.9/22 

M56 

0.4/0. 6 

10/15.2 
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higher output impedance. Dc operating point and common mode rejection are 
obtained by two overlapping feedback loops connected to the current source of 
each differential pair. 

From the outputs of the second stage to the outputs of the first 
stage are two compensation capacitors. These capacitors make the second stage 
look like an integrator, creating a dominant pole in the amplifier's frequency 
response. This effect is identical to the effect of the compensation capacitor 
in )j,A7^1"type op-amps. The open loop gain of the amplifier is then the gm of 
the input differentia] transistors M3 and M4, times the impedance of the compen- 
sation capacitors. There is naturally some attenuation through the two output 
buffers . 


m. 


A = O' 


BUF SC 


COMP 


where 

O' = attenuation of buffer, 

BUr 

q = forward transfer conductance ratio of MOSFET M3, and 
m^ 

CcoMp ^ value of compensation capacitor. 

This formula will hold true for frequencies below the crossover 
frequency of the inner loop formed by the buffers Ml4 and M15, the second gain 
stage, and the feedback path through the compensation capacitors back to the 
buffers' input. Therefore, the crossover frequency of the op-amp with its 
integration capacitors must be well below the unity gain point of the inner 
1 oop. 


The common mode feedback loop dynamics are discussed by stages. In 
the first stage the gain is the gm of Ml (the first stage current source) times 
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the impedance at nodes 7 and 8 (the first stage outputs) in parallel. The 
dominant pole is created by the capacitance at those nodes, including the 
compensation capacitors. The other ends of the compensation capacitors (nodes 
18 and 19) tend to be held at virtual grounds by the action of the second stage 
common node feedback loop. 

^ = "bUF 9m, ^ . 


where 

Qf_|||_ = attenuation of buffer, 
dUF 

g = transconductance of MOSFET Ml , and 

mi 

= compensation capacitor. 

The action of the second stage common mode loop is similar to that 
of the first stage. Complications arise due to common mode feedback through 
the compensation capacitors and interactions with the first stage feedback loop. 
The effect of the first interaction is to create an unwanted pole due to the 
capacitance at node 15, the sources of the second stage differential pair. 

Though the second interaction tends to counteract this effect, it is still 
advisable to minimize the capacitance at node 15 by using the buffer M53 to 
drive Ml , 


As in the first stage feedback loop, care must be taken to prevent 
excessive phase shifting in the buffers which level-translate back down to 
the current source FETs gate. 

The computer-pred icted parameters of the two common mode loops are 
listed on the following page. 
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Low Frequency Gain 
Crossover Frequency 
Phase Margin 


First Stage 
5.5 
14 MHz 

67° 


Second Stage 

5.5 

22 MHz 
60 ° 


Each amplifier requ i red 47 mW dc power and occupies approximately O.658 mm' 
(1020 mil^). 


Both the input amplifiers and the DC I amplifiers were modeled using 
the block diagram of Figure 22 . From Figure 22 the closed-loop transfer 
function is 


V^(s) a (s) 

“ V. (s) ~ 1 + a(s)f(s) 

where a(s) is the open loop gain, 

f(s) is the feedback gain, and 
a(s)f(s) is the loop transmission gain. 

The stability of each amplifier was based on the loop transmission, a(s)f(s). 


A schematic of the input differential amplifier Is shown in Figure 
23. Transistors and Q2 ^^e used to provide the proper input bias level for 
the amplifier. Assuming the open loop gain is infinite, the closed loop gain is 


A(s) 



3.0 pF 

0.5 pF " ^ ’ 


where C,., is the value of the input capacitance and 
I N 

Cpg is the feedback capacitance. 

Both the input capacitors and the feedback capacitors were formed using first 
and second polysilicon as the capacitor plates. The dielectric constant of 
the gate oxide between the two polysilicon levels determine the capacitance. 
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The amplifier closed loop gain is determined by a capacitance ratio that is 
also the ratio of the areas of the two polysilicon capacitors. Since these 
areas can be accurately controlled, the gain is very accurate and predictable. 

A gain of 6 was chosen to compensate for the attenuation of the MDACs. Based 
on computer modeling, 1 pF compensation capacitors were chosen. 

A plot of the predicted open loop characteristics of the input amplifier 
is shown in Figure ?M. From this plot the following parameters can be obtained: 


Crossover Frequency 

6 MHz 

Feedback Ga i n 

1/6 

Gain Bandwidth Product 

36 MHz 

Phase Margin 

63° 


The predicted closed loop frequency response is shown in Figure 25. The -3 dB 
frequency is 9 MHz. 

The design of the DC I amplifier is the same as the input amplifier, 
but instead of providing only signal inversion and gain, the DC I amplifiers 
function as integrators, integrating the CCD output signal onto the feedback 
capacitors. A model of the DC I circuit is shown in Figure 26. The CCD output 
is shown as the split electrodes with a clock driver for each side of the 
electrodes. Two series-pass MOSFETs and Q. 2 ) sre shown that gate the 
signal into the DC I amplifier. The DC! amplifier uses the 0^ pulse as a reset 
pulse. The timing of the DC I circuitry is shown in Figure 27. 

To prevent the DC I amplifier from seeing a large common mode signal 
(CCD 0*)» series-pass transistors are used to isolate the amplifier from the 
split electrodes. To keep the amplifier outputs from drifting apart when the 
amplifier is isolated, M7 and M23 (Figure 2o) were added. When the reset 
pulse goes high, the outputs of each differential stage are shorted together, 
holding the outputs of the amplifier at the same bias point. When the reset 
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Open Loop Ga i n 



Frequency (Hz) 

Figure 24. -D i fferent ia I Amplifier Theoretical Open Loop 




Phose 


(ddjBdQ) aseqd 
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(9p) uieQ 


Frequency (Hz) 

Figure 25.- Input Amplifier Closed Loop Frequency Response. 
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Figure 26.- DCI Amplifie 




















pulse is off and the DCI sample pulse is high, the amplifier outputs settle 
to a new bias point and the signal output of the amplifier is the difference 
between this new bias point and voltage swing of the amplifier. 

+ 

While the CCD (t)^ clock lines are held low, the DCI amplifier is reset 

with the DCI reset pulse (0^). Also, the DCI sample pulse is off, isolating 
+ + 
the 0^ electrodes from the DCI amplifier. Before the 0^ electrodes are allowed 

to float, the DCI sample pulse goes high, connecting the amplifier to the split 

electrodes. Next, the reset pulse goes low, and the amplifier is ready to 

integrate the signal charge onto its feedback capacitors when the CCD 0^ clock 

goes low. A photograph showing the CCD clocks Is shown in Figure 28. The 
+ 

time the CCD 0^ clock is floating is indicated by the arrows. 

The gain from the CCD filter input to the DCI amplifier output can 
be calculated from the input charge of the CCD filter 

Q.. = C. V. 

in I n I n 


where C. is the input well capacity and 
in 

V.^ is the input voltage, 
and the output charge 


a „ = N(i) C. V. 
out in in 


where N is the number of points in the transform (32). 

The factor of arises when a non-dc input is assumed. Since 


Q = C V 
out out out 


then 


out 


^out 


out 


16 C. V. 
in I n 

C ^ 
out 


where C ^ is the feedback capacitor of the DCI amplifier, 
out 
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•- CCD 02 Floating 



Figure 28. -CCD Clocks. 


k8 


■ ■■III 


III 


I I I 



Solving for C 


4.,.pp(^) . 

To provide some gain to overcome the attenuation of the following 
ac-coupl e/dc-restore circuitry and the output source-follower, a value of 
1.4 pF was chosen for the feedback capacitor. Substituting this value into the 
gain equation yields 


out 


16 C. V. 


in I n 


out 


V 


= 16 X 0.26 


ou t 


/V. X 

- 


V 


Gain = 


out 


16 C. 16x0. 26 pF 


I n 


V 


out 


1 .4 pF 


= 3 


A model for computer simulation of the DCI amplifier only is shown 
in Figure 29. This model was constructed to examine the characteristics of 
the amplifier. The CCD split electrode has been modeled as an equivalent 
capacitor, and the on- res i s tance of the series-pass transistors has been cal- 
culated to be 160 0, based on the physical size of the MOSFETs and the gate-to- 

+ 

source voltage. A model for calculating the CCD 0^ capacitance is shown in 
Figure 30. The values of the capacitors for two filters in parallel are: 


Depletion capacitance 

Oxide capacitance 

First poly-to-moat and 
metal -to- substrate 

Overlap of CCD 0^ 

Overlap to CCD0^ 

These capacitances can be combined Into one 
ca 1 cu 1 a t ions . 


C^ = 4.6 pF 
dep 

Cox 

= 5.38 pF 


C^ = 19.2 pF 
C^ = 19.2 pF . 
equivalent capacitor for gain 
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Figure 29.-Model for DCI. 



The predicted open loop characteristics of the DC I amplifier are shown 
in Figure 31- From this figure the following parameters can be calculated: 


Crossover frequency 17 MHz 

Feedback gain 1/32 

Gain-bandwidth product 550 MHz 

Phase margin 65°. 


The predicted closed loop frequency response of the DC i amplifiers is shown in 
Figure 32, and the -3 dB frequency is seen to be 25 MHz, 

2 . Experimental Results 

Experimental results have been obtained for the input amplifiers 
using the test setup shown in Figure 33- A pulse generator was used as the 
inputs to the amplifiers, and attenuators were used to obtain low-level 
signals. Because there are two input amplifiers, the inverting inputs 
and the noninverting inputs were connected to allow testing of both ampli- 
fiers. The input signal and the inputs of the amplifiers were connected 
through the connection box in a number of ways, depending on the test performed. 
Because the output signal of the amplifiers has a large dc component, an ac- 
coupl e/dc-restore circuit was used. Finally, the output signal was sampled so 
that a digital voltmeter could be used to record the test results. As can be 
seen in Figure 33» both the input signals and all the output signals were 
recorded through the same ac-couple/dc-restore and sampl e-and-hold circuit. 

The results of the amplifier tests are shown in Figures 3^- - 36. 

Figure 3^- shows the single-ended transfer characteristics of the amplifier. 

The single-ended gain is approximately 3 , and output is linear up to about 
0.4 V input. Figure 35 's a plot of the response of the amplifier to a 
common mode input signal. The differential mode characteristics are shown in 
Figure 36 - 


52 








I 














ristics of Input Differential Amplifier. 



Figure 35.- Response of Input Differential Amplifier to Common Mode Input 








Some photographs of the input amplifier in operation are shown in 
Figures 37 and 38 . Figure 37 shows the single-ended output of the amplifier 
with a ramp input. Figure 38 shows two photographs of the input amplifier. 

In the top photograph the input and output signals of the amplifier are shown, 
and again the gain can be measured to be approximately 3. The bottom photo- 
graph is better for determining gain and linearity. The horizontal axis has 
been calibrated so that it corresponds to the differential Input voltage. The 
vertical axis is the single-ended output of the amplifier. By comparing the 
amplifier output to the ideal straight-line output which has been superimposed, 
one can see the output becomes very nonlinear above 0.6 V input. The open-loop 
transient characteristics of the DC I amplifier are shown in Figure 39. The 
top photograph shows the open-loop dc transfer curve, and the bottom photo- 
graph shows the open-loop step response. The rise time is about 1.6 ps . 

D . Multiplying D i g i ta 1 -to-Ana log Converters (hPACs) 

Four multiplying D/A converters (MDACs) are on the CZT 1C. These perform 
the pre- and post-multipliers of the CZT algorithm. They also perform multi- 
plication by the reference for the correlation and filtering configurations. 
Each MDAC output is an analog signal whose value is proportional to the product 
of an analog input signal and a digital word. For the CZT pre- and post- 
multiplications, the required chirp signals are stored digitally on chip in a 
ROM. For multiplication by a reference, the digital signal must be supplied 
externa 1 ly . 

1 . Des i qn 

Each MDAC consists of an input sample-and-hold circuit, ROM sense 
and MDAC drive circuitry, a capacitor array and an output buffer. A photo- 
graph of two MDACs is shown in Figure 40. 

The MDAC minus the sample-and-hold and ROM sense/MDAC circuitry is 
shown in Figure 4l with W/L and capacitor values indicated. This circuit 
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out 



Both Traces 
200 mV/d i V 
Input and Output 


I nput 


3.0 

2.5 

2.0 

1.5 
1 . 0 
0.5 

0 



Ideal Output 

Output of Input 
Ampl i f i er 


0 


0.4 0.8 

V. (V) 
in ' ' 


Figure 38 . - Photograph Showing the Gain and Linearity 
of Input Amplifiers. 
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(a) 



DC I Ampl i f ier 

Output 

0.5 V/div 


Input (XlOO) 
0.5 V/div 
1 jj.s/div 


(b) 


Figure 39.-DCI Amplifier Open Loop Cha racter i s t i cSc 
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performs a full four-quadrant multiplication using a capacitor divider network. 

For noninverting operation, a sampled input signal is applied to the V line, 

and a reference voltage is applied to the V line. During the initial portion 

of each clock cycle, the common output node of the capacitor array is 

clamped to through MI50 by (refer to Figure 44), Also, the switching 

signals v. are high, and v. are low. (The v. and v. waveforms are nonoverlapping 

complementary waveforms.) This connects the largest capacitor to the signal 

and all other capacitors to the reference V through the switching transistors 

M205 through M275. Next, 0^ goes low, leaving the common node floating, A 

complementary pulse to (0'p is capacitively coupled to the common node by 

MI5I to reduce the level translation of V ^ when MI50 Is turned off. Next, 

ou t 

.the V. and v. signals correspond i ng to "l"s at the i bits, b., of the 

digital word are inverted, causing capacitors C. to be switched from one input 

line to the other through the M2i6s. The v. and v, corresponding to B. = "0" 

remain unchanged. The result of switching capacitor Cj from one Input line 

to the other is that the common node voltage, V , is changed by an amount 
j ^ out 3 / 

± (V - V ) r (C /2 X b,)/(2C + C )]. This change has a positive sign for 

' + o I o stray 

1=1 to 7 and a negative sign for i = 0, since b^ is the sign bit. The 
resulting output, for multiplication by an arbitrary digital word is 

given by superposition to be 


V 

ou 


t 


V . +(V 
ref + 


-vj( 


2C + C 


stray 


n 

[-b + S b. 2"'] . 

° i = l ' 


Two's complement coding is used for the digital word. 


The output voltage of the MDAC is equal to the product of the 

signal voltage (V_^ “ V ) and the digital input word times the attenuation 

factor C / (2C + C ). For C =0, the maximum possible gain is i. 

o ^ o stray^ stray t' ^ ^ 

The actual stray capacitance was calculated to be 5 pF, making the gain to the 

common node V = 10.56/(2 x 10.56 + 5) = 0.4. The capacitor array was formed 
out 
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by using first polysilicon to second polysilicon capacitance. First poly- 
silicon was used as the output node, and the major source of the stray capaci- 
tance is from the first polysilicon to the substrate. To control the capaci- 
tance ratios in the capacitor divider network, the larger capacitors were 
formed as binary multiples of a smaller capacitor. Capacitor is a 
51 pm X 5 1 pm (2 mil x 2 mi 1 ) poly-poly capacitor with a value of about 0.66 pF 
Capacitor C^ Is obtained by tying two 51 pm x 5 1 pm poly-poly capacitors 
together so that C^ = 2C^^ = 1.32 pF. The three larger capacitors are formed 
by tying 51 pm x 51 pm capacitors together (C^ = = 8C^, C^ = 16C4), 

In this way, the correct capacitor ratios are obtained without maintaining 
perfect control of the absolute capacitance. The capacitors C^ , C^ , and C^ 
are formed by smaller poly-poly areas, since a small error in their value will 
not significantly affect the output accuracy. An attempt was made in layout to 
keep all 51 pm x 5 1 pm capacitors identical by placing dummy first poly around 

certain parts of the array to maintain identical first poly edges. The V and 

V inputs required for the MDAC are supplied by the buffered samp 1 e-and-hol d 
circuit shown in Figure 42. The v"*" and V outputs are from identical stages 

so that the circuit can be used in either an inverting or a noninverting mode. 

For noninverting operation, the signal voltage is applied to the + input, and 
a reference voltage corresponding to the dc value of the signal voltage is 
applied to the - input. The signal is buffered by the first souce-fol 1 ower 
(MlOl, M102, M103) and is sampled by the input samp 1 e-and-hol d pulse through 
M104. The holding capacitor consists of a 2 pF first polysilicon to second 
polysilicon capacitor and the gate capacitance of MlO/ for a total of about 
5.3 pF. The sampled signal is buffered by a second source-follower (M105, M106 
MIO 7 ) before being applied to the MDAC. This second source-follower was 
designed to drive the capacitive load of the MDAC at clock frequencies of up 
to 2.5 MHz. This accounts for the large device size of MIO 7 . 0^ is a pulse 
approximately complementary to the input samp 1 e-and-hol d pulse. It is applied 
to the gated capacitor Ml 14 to reduce the nonlinear level translation occurring 
when M104 is turned off. Dual series current sources are used In the source- 
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followers to decrease the effect of A. (the channel length modulation parameter) 
on the gain. The gates of the current sources are fed by the voltage division 
stage consisting of Ml08, M109f and Ml 10. The gain of each source-follower 
is 0.9, making the gain of the entire buffered sampl e-and-hol d circuit 0.81. 


The switching signals v. and v. are supplied by the ROM sense and 

MDAC drive circuitry shown in Figure 43- Each MDAC contains eight of these 

circuits. , 0 and 0 are complementary nonoverlapping clocks. 0 comes up 
Be B 

during the first portion of each clock cycle, pulling v. up (if it is not already 
up), and charging the gates of M201 and M204 to the value of the ROM (or exter- 
nal word) output voltage. This voltage will be high if b. = "1" and low if 
b, = "0". If b. = "1", the gate of M201 will be high and v. will fall with 
0 . Then 0 comes up, pulling v. up since the gate of M204 is high. Ml 14 is 
a gated capacitor that bootstraps the gate voltage of M204 high enough to 
allow V. to reach the full voltage of 0^. At the end of the clock cycle, v. 
falls with 0^. If b, = "0", the gate of M201 will be low, and v. remains high 

when 0 falls. Since the gate of M204 is uncharged, v, does not come up with 
B I 

<*c- 


The switching waveforms and typical signal waveforms for a multi- 
plication by 0.5 are illustrated in Figure 44. The digital word is 01000000. 

The Vj and v, pulses switch the bottom plate of capacitor C, from the V line to 
the V line (Figure 41) to accomplish the multiplication by 0.5- The MDAC out- 
put is valid while 0^ is up after the signal has settled. 

The only dc power consumed by the entire MDAC is in the four source- 
followers and the voltage divider of the buffered sampl e-and-hol d stage and 
in the output buffer of the MDAC. The total dc power consumed by these stages 
in one MDAC is 40 mW. 


The total attenuation of the signal from the buffered input to the 
buffered output is 0.3. This is the result of transmission through three 
source-followers with gains of 0.9 each, and the MDAC capacitor array with a 
gain of 0.4. 
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2 . 


Experimental Results 


Figure 45 is a plot of the transfer curve for each side of the input 
sample-and-hold of Figure 42. The gain is 0.84, and the nonlinearity is very low 
(< 0.2% of input span). The input sample-and-hold circuitry functions properly 
for inputs of 1 to 7 V. 

The MDACs were evaluated using the test setup shown in Figure 46. 

The MDAC was tested with both an external digital signal and chirps stored in 
internal ROM. 

The evaluation of the MDAC using the sine chirp of the ROM is shown 
in Figure 47 and Table II. Figure 47 shows the response of the MDAC to a dc 
input and the sine chirp. Table II shows the word number, its binary code, 
and the analog value of the binary code for the sine chirp. The MDAC output 
voltages were normalized and compared to the analog representation of the code 
stored in the ROM. 

The errors were found to be not more than ± 0.25 least significant 
bit (LSB), or one part in 512. Since the ROM is accurate to within ± 0.5 LSB, 
the worst-case error from the MDAC should be ± 0.75 LSB. 

The MDAC was also tested using external digital words, and the 
relative accuracy is plotted in Figure 48. The differential nonlinearity was 
measured to be 0.25 LSB. Other parameters measured on the MDACs using external 
digital words were digital feedthrough , analog feedthrough, analog nonlinearity, 
and analog differential nonlinearity. The results of these tests are summarized 
i n Tab 1 e III . 

E. Read-Only Memory (ROM) and ROM Address 
1 . Read-Only Memory 

The on-chip ROM that supplies the cosine and sine chirp signals to 
the MDACs is illustrated in Figure 4^. Two ROMs like the one shown are on the 
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Figure 45. -Transfer Curve for Input Samp I e-and-Ho 1 d 
C i rcu i t of MDAC. 
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Figure 47. -Response of MDAC to dc Input and Sine Chirp of Internal Rom 
(V+ - V) < 0. 
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TABLE II 


MDAC ERROR 


Word 

Number 

B Inary 
Code 

Analog 
Va 1 ue 

MDAC 

Output'' 

Error‘s 

0 

00000000 

0.0000 

0.0000 

0 

1 

00001 100 

0 . 0945 

0.0951 

0.077 

2 

oonoool 

0.3858 

0.3870 

0. 154 

3 

01100010 

0.7717 

0.7725 

0.103 

4 

on n n 1 

1.0000 

1.0000 

0.000 

5 

01010001 

0.6378 

0.6388 

0. 128 

6 

11001111 

-0.3858 

-0.3843 

0.192 

7 

10000010 

-0.9921 

-0.9933 

0.154 

8 

00000000 

0 

0.0009 

0. 1 15 

9 

01111110 

0.9921 

0.9920 

0.013 

10 

nooiin 

-0.3858 

-0.3844 

0.179 

1 1 

lOlOl 1 1 1 

-0.6378 

-0.6373 

0.064 

12 

on 1 n 1 1 

1 . 0000 

1 . 0004 

0.051 

13 

1001 1 1 10 

-0.7717 

-0.7715 

0.026 

14 

001 10001 

0.3858 

0.3872 

0. 179 

15 

11110100 

-0.0945 

-0.0938 

0. 090 

16 

00000000 

0,0000 

-0.0002 

0.026 


Norma 1 i zed 

+ 

Fraction of LSB 
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-LSB 
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Figure 48 .- Relative Accuracy of MDACs. 



TABLE III 


MDAC CHARACTER 1ST ICS 


Analog Nonl inea r i ty ^ 

( 2 ) 

Analog Feedthrough 

Analog Differential Nonlinearity 

(4) 

Digital Feedthrough 


(3) 


0.2% of input span 
< 0.5 niV at 100 kHz 
± 0.25 LSB 
± 1.25 mV 


( 1 ) 

( 2 ) 

(3) 

(^) 


Measured by holding digital input constant and varying the analog 
Variable analog input with digital word programmed for zero gain. 
Maximum deviation of any bit size from theoretical value of 1 LSB 
Variable digital word with analog input held at zero. 


i npu t . 
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1C. One supplies the cos chirp and the other supplies the sine chirp. The 
outputs b. are attached to the ROM sense circuit of Figure 43- These outputs 
can come from an external digital input by holding the external select line 
high and the internal select line low. To use the internal ROM, the internal 
select line is held high and the external select line is held low (to reduce 
capacitive loading of the ROM). There is actually a second set of internal/ 

external select switches and outputs b. on chip. Each set feeds a pair of MDACs. 

' 2 
A photograph of the ROM is shown in Figure 50. It occupies an area of 0.62 mm 

(35.5 X 27 mil^). 

The ROM shown consists of 17 eight-bit words. One word consists of 
a horizontal row of eight possible transistors in Figure 49. A bit is coded as 
a "0" or a "1" by the presence or absence, respectively, of a transistor. The 
ROM address lines a^ are fed by a 17-stage up/down ring counter. The ROM is 
addressed In an up/down manner since the 32-point chirp is symmetrical about 
the seventeenth point. Each word, j, is addressed sequentially by pulling the 
appropriate address line to a high voltage. All other address lines are held 
low. pulls the output lines b. high, and when 0^ falls, the output b. 

will fall if there is a transistor in bit i of word j. if there is no tran- 

sistor there, b. will remain high. For example. If word two (line a^) is 
addressed in Figure 49, outputs b^ and b^ will be low "0," and outputs b^ and 
b^ will be high "1". To illustrate the ROM timing, refer to Figure 44. 

The appropriate address line, a^ , comes up and down with 0^^ so that 
aj is high when 0^ falls. The ROM output is valid only after 0^ has fallen. 

The fact that all ROM outputs are high when 0^ is high insures that M204 of 
the ROM sense circuit will discharge the gate of M206 of Figure 41 once each 
cycle so that no excess charge can build up on this node. 

A plot of the sine chirp coded in the ROM is shown in Figure 51. 

This plot can be compared to the photograph in Figure 47. 
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Figure 51-'Analog Representation of Sine 




2 . 


Up/Down Ring Counter (ROM Address) 


The ROM address lines a^ - a^^ are fed by a l7"Stage shift register 
that has the ability to reverse the shifting direction. A three-stage version 
of this circuitry is shown in Figure 52. Each basic shift register stage (one 
bit) is outlined in Figure 52. 0^ and 0^ are the nonoverlapping two phase 

clocks input to the 1C. The basic shift register consists of two dynamic 
inverters with transmission gates, and its operation depends on charge storage. 
The operation of shift register 0 (SRO) of Figure 52 is illustrated in Figure 
53 . When 0 ^ comes up, the 0.2 pF capacitor on node 2101 is charged to the 
level of the input to SRO through M2100. Also, the 1 pF capacitor on the first 
inverter output node 2102 is charged to a high voltage through M2 102. When 
0j falls, node 2102 will discharge through M2101 if node 2101 is high; other- 
wise, node 2102 remains high. Thus, node 2102 is the inverse of the input to 
SRO. Next, 0^ comes up, turning on M2 1 03 . If node 2102 is high, the charge on 
the 1 pF capacitor will divide between itself and the 0.2 pF capacitor on node 

2103 , and the voltage on node 2103 will be [1 pF/(l pF + 0.2 pF)] v^, where 
Vq was the voltage on node 2102 before 0^ turned M2 103 on. If node 2102 is 
low when 0^ comes up, any charge on node 2103 will discharge through M2 101, 
since the gate of this transistor must be high. 0^ also charges the 1 pF 
output capacitor on the second inverter output node 2104 through M2105. When 
02 falls, node 2104 will remain high if node 2103 is low or will discharge 
through M2 1 04 if node 2103 is high. Node 2103 is now isolated from 2102, since 
M2103 is off, allowing node 2102 to change without affecting the output node 

2104. Thus, the voltage at the output of SRO is that of the input delayed by 
one clock cycle. 

To address the ROM properly, the shift register must shift one "1" 
back and forth across its 17 stages. This requires circuitry to detect the 
" 1 " when it reaches each end of the shift register and to reverse the transfer 
direction. It also requires circuitry to insure that only one stage of the 
shift register contain a " 1 ," and that the other sixteen stages contain " 0 "s. 
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Figure 53.-Typicai Shift Register Operation (SRO) . Capacitive coupling, time delays, 
and MOSFET thresholds ignored. Previous SRO inputs have been "0"s. 
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The shift register must also have circuitry to synchronize it with an external 
pulse. 


The direction of transfer is determined by the state of the static 
S-R flip-flop FFl. With flip-flop output node 2182 high and node 2183 low, the 
feed-forward transmission gates M2107, M2106, and M2116 are on, and the feedback 
transmission gates M2127j M2126, and M2117 are off. The output of each stage 
of the shift register is electrically connected to the input of the next stage, 
and transfer is in the forward direction. When the "1" reaches the final stage 
of the shift register, transistor M2128 will be on and the reset node 2180 
will come up and down with 0^ This causes the flip-flop to change state, 
and the feed-forward transistors are turned off while the feedback transistors 
are turned on. Each stage's output is now electrically connected to the pre- 
vious stage's input. Node 2111 rises to a high voltage, since it is now 
connected to node 2124. The "1" appears next at the output of the next-to-last 
stage (SR 1 in this example). "0"s are loaded into the last stage by M2127, 
which is tied to ground. Transfer continues in the reverse direction until 
the "1" reaches the output of the first stage. Then when 0^ comes up, M2 1 08 is 
on and the set node 2181 comes up and down, causing the flip-flop to change 
stage again. The feed-forward transistors are now on, and the feedback tran- 
sistors are off. The "1" next appears at the output of the second stage. "0"s 
are loaded into the first stage by M2107, which is connected to ground. 

To insure that the shift register contains only one "1," two cases 
must be considered. The first case is that when power is supplied to the 
circuit, the shift register contains one or more "l"s. In this case, the out- 
put of the first stage of the shift register will eventually be high. When 
this occurs, M2 1 08 will be on, and node 2181 will come up and down with 0j . 

This sets the inputs to each stage of the shift register low through M2109 
and M2 129, except for the second stage, which is set high by M2 11 9. At the 
end of this clock cycle only the second shift register contains a "1," and 
transfer is in the forward direction. 


85 



The second case is that when power is supplied to the circuit, the 
shift register contains all “0"s. In this case, a 17“input dynamic NOR cir- 
cuit sets the first stage to "1." The NOR circuit consists of transistors 
M2130 - M2133, in Figure $2, Node 2130 is pulled high with 0^ If the input 
to any stage is high, node 2130 falls with 0^, since at least one of the tran- 
sistors M2130 - M2132 Is on. If the inputs to all stages are low, then all 
transistors M2130 - M2132 are off, and node 2130 will remain high. This leaves 
M2 134 on. Node 2131 will then come up and down with 02* This discharges node 
2103 through M2135 so that when 0^ falls, the output of the first stage, node 
2104, will remain high or "1." After this occurs, the shift register contains 

a "1," and node 2131 should never come up again. Node 2130, however, will 

rise and fall with 0^ every clock cycle. 

Three transistors, M2138, M2139, and M2191, are added to the circuit 
of Figure 52 to synchronize the counter with a sync pulse (M2 1 1 9 is also only 
necessary for synchron ization) . The sync pulse is derived from another circuit. 
When it occurs, it comes up with 0^ and down with the rising edge of cf)^. The 
set node 2l8l is pulled up with 0^ since M2138 is on. This resets the inputs 
of all stages except the second low through M2109 and M2129. The second stage 
input is set high by M2119. Also, node 2103 is discharged by M2139 so that 
M2 1 04 will not sink current from node 2111 to ground. The sync pulse is also 
applied to the gate of M2191 to make sure that the flip-flop comes up in the 

forward transfer state (node 2182 high). This transistor is necessary because 

it is possible for the reset pulse to occur when a "!•' is in the last stage, 
in which case the reset line 2180 would also come up-and down with 0^ Since 
the sync pulse remains high after 0^ is low again, the flip-flop will be set 
in the proper stage. At the end of the clock cycle, the output of the second 
stage will be “1," and the shift register will be transferring in the forward 
direction. (M2191 is necessary only to guarantee synchronization after the 
first sync pulse. This is not required here. If M2191 were removed, synchroni- 
zation would be guaranteed after the second sync pulse.) 
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The shift register outputs go to the gates of transistors M2150 - 
M2 152 , which allow 0^ to pull up the ROM address line attached to the transis- 
tor whose gate is high. The depletion devices (M2160 - M2162) are attached to 
the ROM address lines to keep capacitive coupling in the ROM from pulling the 
other address lines up. The depletion device on the address line that is high 
will sink 0,6 mA of current from 0^^ . The 0^^ driver is able to supply this 
current with very little drop in output voltage. The output pulse is illus- 
trated in Figure 53. Figure 5^ illustrates the counter in normal operation 
(i.e., one "1" bit, counter synchronized to sync pulse 3004 if it exists). 

It shows the first and last shift register outputs, the first and last ROM 
outputs, the relationship to the sync pulse, and the S-R flip-flop Inputs and 
outputs. The flip-flop FFl consists of a pair of static inverters with boot- 
strapped loads. M 2 I 87 and M2188 keep the gates of M2 1 85 and M2186 from falling 
below a threshold below M2 1 89 and M2 1 90 are gated capacitors that boot- 

strap the gates of M2185 or M2186 when the output 2182 or 2183 starts to come 
up. This turns the load devices on hard, decreasing the rise time and allow- 
ing the output voltage to reach The set and reset transistors M2183 and 

M2 182 are in parallel with the driver devices M2 184 and M2l8l. 

A photograph of the up/down counter is shown in Figure 55. Its 

2 2 

area is approximately 0.75 mm (40 x 29 mi 1 ). 

3 . Synch r on i za t i on 

To operate two or more ICs together, the up/down counters must be 
synchronized to each other. it is required that the counter of an 1C be 180° 
out of phase with the counter it is synchronized to. One cycle of the chirp 
z-transform algorithm consists of 64 clock cycles. Therefore, when two ICs 
are operated together, one should be on count 32 when the other is on 

count 0 (t^). 
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(tj,) (t32) (t33) 



Figure 5^. “Up/Down Counter Waveforms - Norma] Operation (Sychronized to 3004) 
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Figure 55. •” Photograph of Up/Down Counter. 



In order to obtain this synchronization, each counter must have an 
external sync pulse output to provide a pulse to the other IC, and a sync 
pulse input to accept a pulse for synchronization to another IC. The output 
sync pulse is a square wave. The sync pulse input c i rcu i t accepts the square 
wave output from another IC and converts this to an internal pulse (3004) 
with duration of less than one clock cycle to reset the counter and external 
sync pulse generation (refer to Figures 52 and 54). The relationship of the 
sync pulses for two synchronized chips is illustrated in Figure 56 - Also 
shown is the relationship of the sync pulses to the flip-flop FFl (Figures 52 
and 54 ) output (2182) to the feed-forward transistors in the up/down counter. 

It can be seen that the up/down counters of the two I Cs that have a period of 
32 clock cycles are in phase, and the sync pulse output of the second IC is 
180 ° out of phase with that of the first stage. The external sync pulse 
generator is shown in Figure 57. The circuit uses a set/rest flip-flop 
(identical to FFl in the up/down counter) and seven more transistors to imple- 
ment a toggle flip-flop that changes state each time a pulse appears on the 
set line 2181 to FFl. This occurs once every 32 clock cycles in normal 
operation. The basic RS flip-flop consists of the unmarked transistors in 
Figure 57. Figure 58 illustrates the circuit's operation. Assume the output 
2282 is low and 2283 is high. Then node 2286 will be low and node 2287 will 
be high, since 0^ will have turned M2289 and M2290 on and off. Node 2181 will 
come up and down with 0^ every 32 clock cycles. With the flip-flop in the 
above state, M 2292 is on, and node 2281 will come up and down with the input 
node 2181 . M 229 I is off and node 2280 stays low. This causes the flip-flop 

to change state. The transistor M2293 is on during this process to keep node 
2280 from being pulled up by capacitive coupling to the flip-flop. M2294 is off. 
Now output 2282 is high and 2283 is low. When 0^ comes up, node 2287 dis- 
charges through M2290, and node 2286 charges through M2289- The circuit is 
triggered during synchronization so that the output waveform at 2282 
is 180 ° out of phase with the sync pulse. The duration of the internal 
sync pulse (3004) is longer than that of 2181, so that sync is assured. The 
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Figure 56.-Synchron izat ion of Two CCD CZT ICs* 
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relationship of the internal sync pulse (3004) to the generated external sync 
pulse is shown in Figure 58 during normal operation. It is out of phase with 
node 2282 . The actual sync pulse generated for use off-chip comes from a 
driver that will be discussed later. The inputs to the driver are 2282 and 
2283 , and the output waveform is that of node 2282 with slightly more rise and 
fall time delays. This driver is necessary to drive the larger capacitance 
associated with going off-chip. The output pulse is called T^ • 

The circuit used to generate the short internal sync pulse from an 
applied external sync squarewave is called GEN30 and is shown in Figure 59« 
Figure 6 O illustrates its operation. When the external sync input is low, 
node 3001 remains high after falls. Node 3002 comes up and down with 0^ 
every clock cycle, and node 3003 is charged to a threshold below When 

the external input comes up (which will be with 0^ plus a small delay), the 
output node 3004 comes up immediately, since M3003 is on. The gate of M3005 is 
bootstrapped by M3008 so that 3004 reaches the full potential of the input. 

Node 3001 falls immediately with the input rising edge, leaving 3002 floating 
high. When comes up again, node 3001 comes up, allowing 3002 to discharge 
through M3003. M3006 and M3007 are turned on, turning off M3005 and pulling 

the output node 3004 down. The output cannot come up again until M3004 charges 
node 3003 high. This cannot occur until after the input external pulse falls. 
In this way, GEN30 generates a pulse with a duration of about g clock cycle 
each time the external input comes up. 

Figure 6l is a photograph of the sync outputs of two CZT ICs 

operating together. The sync output of one 1C is fed into the sync input of 

* 

the other. The sync output waveforms have a period 64 times that of the 10 kHz 
input clocks, 0^ and (i)^. Note that the two sync outputs are 180° out of phase 
as in Figure 58- 
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Figure 60 -GEN 30 Timing^ 



Figure 6 1.- Synch ron i za t ion Outputs of Two CZT ICs 
Operating Together- 


F . Timing Circuitry 

1 . Clock Generators 

The derivation of the internal clocking pulses from the two-phase 
inputs 0| and 0^ is illustrated in Figure 62. The 13 generators shown there 
use one or more of the circuits described below. 

The first circuit will be referred to as Gl and is illustrated in 
Figure 63. This circuit configuration is used to produce most of the clock 
pulses. It requires a pull up pulse (node 4000) and a nonoverlapping pulldown 
pulse (node 4001). Assume that a pulldown pulse has occurred and the output 
(4006) and all internal nodes are low except for node 4003, which is high. 

When the pullup pulse occurs at node 4000, nodes 4002 and 4004 closely track 
the Input. Nodes 4005 and 4006 are slightly delayed because 4003 is still 
high, keeping m 4006 and m 4008 on. Since 4005 is held down and 4004 rises 
quickly, the gated capacitor M4010 is charged to a level of - V^. When 
4002 has risen high enough to turn M4002 on, pulling node 4003 down, m 4006 
and M4008 are turned off. At this point nodes 4005 and 4006 begin rising 
rapidly, and node 4004 is bootstrapped by m 4010 to a level 2 This 

increases the rise times and allows the output node 4006 to reach The 

time delays for the actual circuits used are typically 20 ns. The risetimes 
are typically 20 ns to 30 ns. When the pulldown pulse occurs at node 4001, 
node 4002 is pulled down, allowing 4003 to be pulled up while node 4004 is 
pulled down. Nodes 4005 and 4006 are then pulled down. The time delay before 
node 4006 starts down is typically 10 ns, and the fall time is typically 10 ns 
to 15 ns. The device sizes required for the Gl generator to produce 20 ns 
risetimes depend on the capacitive load it must drive. Typical designs set 
the W of m 4009 in mils equal to the capacitive load in pF and the W of m 4008 
was set to 1.5 times that of M4009. The values of the other transistors depend 
on internal node capacitance ratios. Therefore, each Gl driver required a 
custom des ign. 
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Figure 62. -Block Diagram of Clock Generators. 






















A special version of the Gl driver is required to allow the CCD 
electrodes to float. This driver, G2 , is shown in Figure 64. The difference 
between G2 and Gl is that the gate of m 4004 on G2 is tied to a third pulse, 
which comes up when it is desired to float the outputs. The G2 driver also 
has a second pair of output transistors so that one generator drives two of 
the four CCD outputs. A second G2 generator drives the other two. The 
line to the second G2 generator can be trimmed to adjust for offsets in the 
DC is. 


The generator G3 shown in Figure 65 is a bootstrapped inverter circuit. 
The basic inverter consists of M4101 and M4102. m 4103 keeps the gate of m4102 

from falling below a threshold below m 4104 is a gated capacitor which 

bootstraps the gate of M4102 to a high voltage when the output node 4102 comes 
up. This speeds up the risetime by turning m 4102 on hard. The generator G4 
is an inverter whose pullup transistor is not bootstrapped. It is shown in 
Figure 66. 


Figure 67 shows the generator G5 • G5 produces a short pulse at node 
4204 from a longer trigger pulse at node 4200. Before the trigger pulse occurs, 
node 4002 is pulled up by a set pulse on node 4201. When 4200 comes up, the 
output 4204 comes up immediately, since m 4200 is on. The gate of M4200 is 
bootstrapped by gate-drain and gate-source capacitance so that 4204 reaches the 
same potential as 4200. While the trigger pulse is up, a pulldown pulse occurs 
on node 4203, turning off m 4200 and pulling the output down. 


Figure 68 shows 
down pulse are required. 


the simple G6 generator. A pullup pulse and a pull- 
The output is designed to have a maximum voltage of 
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of the generators Gl - G6 used to produce the internal 
Figure 69 . The approximate capacitive load of each 
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Figure 68. -G6 Generator. 
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Figure 69 (b),- Connect ions of the Clock Generators Showing Capacitive Load of 
Each Driver. 
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driver as calculated from the bar layout is indicated. The node numbers corres- 
pond to the names given each set of generators and to the SPICE modeling per- 
formed. Many clock waveforms required series connections of G1 drivers to 
generate propagation delays critical to the chip timing. Figure 70 shows the 
resulting internal timing pulses as predicted by the computer simulation pro- 
gram, SPICE, for a clock period of 400 ns (2.5 MHz operation). 


2 . Experimental Results 

Photographs of some of the actual clock waveforms are shown in 

Figure 71 through 73- These pulses were monitored on-chip with an external 

probe. Figure 71 shows the leading edge of 0 and its pullup pulse 0,. The 

B I 

0 leading edge is delayed from the 0 leading edge through two Gl drivers, as 
B 1 

seen in Figure 69 (a). The actual delay as seen in Figure 71(a) is about 60 ns. 
The risetime is 20 ns. The delay predicted by computer simulation was 55 ns, 
and the predicted risetime was 20 ns, as shown in Figure 70. Figure 71(b) shows 
the trailing edge of 0 and its pulldown pulse 0 . No delay is added to the 
0 trailing edge. The actual trailing edge delay is 30 ns, and the fall time 
is 15 ns. The predicted delay was 10 ns, and the predicted fall time was 15 ns, 
as shown in Figure 70. These results show close agreement between the com- 
puter simulations and the actual generator outputs. The 0 and 0 waveforms 

B L 

were designed to be comp 1 emen tary and nonoverlapping. Figure 72 illustrates 

the success of the design of these clocks. 0^ is generated using the same 

dual Gl driver configuration as 0 with the pullup and pulldown pulses inter- 

B 

changed. Note that-0 and 0 are never high at the same time, and their lead- 

B C 

ing and trailing edges do not overlap each other. Figure 72 may be compared to 
the computer-simulated timing of Figure 70. 


A photograph of two of the larger clock generators (CCD 0~) is shown 
in Figure 73- These are G2 drivers as in Figure 64. It measures 1.75 x 0.4 mm' 
(69 X 16 m i 1 ^ ) . 
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Figure 70. -System Timing. 
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G . Experimental Results 

Evaluation of the CZT 1C has been confined to obtaining power density 
spectra of real inputs using the system shown in Figure ~Jh ■ This system requires 
only two MOS level clocks, voltage bias for the various circuits, and the 
squaring function. Figures 75 through 79 are photographs showing the operation 
of the CZT 1C. Figure 75 shows the operation of the 1C with 10 kHz clocks 
and no postmultiplication. The upper photograph shows the response of the 
MDACs to a dc input and the internal ROM. The bottom photograph shows the 
response of the 1C to a 1.8 kHz sinusoid. The output shows the envelope of the 
real and imaginary outputs, since the input is not synchronized to the chip 
clocks. Similar operation of the 1C is shown in Figure 76 with 1 MHz clocks. 

The input signal is a 440 kHz sinusoid with a dc offset. Figure 76 also illus- 
trates a problem with the tapweights of the filters. The CCD tapweights were 
erroneously shifted one bit, causing the 32nd bit to appear as the first bit, 
while the first bit shifts to the second bit, and so on. This will affect 
only spectral power density analysis in that the output sequence is rotated. 

It will also affect function implementations requiring postmultiplications. 

To perform postmultiplications, a second chip will be required whose timing 
is delayed by one clock period. This error can easily be corrected with photo- 
mask changes. 

Figure 77 shows the power density spectra for three different sinusoidal 

inputs with a dc offset. The output is observed to be delayed as the input 

frequency increases from 625 Hz to 4.0 kHz. Similar operation is shown in 

Figure 81 with 1 MHz clocks and the input varying from 125 kHz to 375 kHz. 

The apparent sidelobes are caused by the limited bandwidth and inaccuracies of 

the squaring devices. The response of the CCD CZT 1C to a 312.5 Hz and a 

625 Hz square wave input is shown in Figure 79. The output can be seen to 
2 

decrease as 1/n (where n is the harmonic index). 














(a) Response of MDACs to dc Input and Internal ROM 



(b) Response of IC to 1.8 kHz Sinusoid 


1 16 


Figure 75.-Operation of CCD CZT IC with 10 kHz Clocks. 
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(a) Response of MDACs to dc Input and Internal ROM. 



(b) Response of IC to 440 kHz Sinusoid. 


Figure 76.-0perat ion of CCD CZT IC with 1 MHz clocks. 
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Figure 77.- Power Density Spectra for Three Sine Waves Obtained Using the 
32-Point CCD CZT at a 10 kHz Sample Rate. 
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Figure 78 -Power Density Spectra for Three Sine Waves Obtained Using the 
32 Point CCD CZT at 1 MHz Sample Rate. 
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Input 625 Hz 
Square-Wave 


Output 


Figure 79 -Response of CCD CZT 1 C to Square-Wave Inputs (10 kHz Clock Frequency). 
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H. Appi I cat ions 

The most straightforward application of the 32-point CCD CZT 1C is to 
obtain power density spectra using the block diagram of Figure ~jh. 

With a few external components, a complex input/output CZT can be per- 
formed with one chip as shown in Figure 80. The external components are 
necessary to disable the inputs when a postmultiply operation is needed for 
the true CZT. Because of this blanking, the system in Figure 80 has only a 50% 
input/output duty cycle. To obtain a 100% input/output duty cycle CZT, two 
chips are needed as shown in Figure 81. This configuration does not require 
the external circuitry of Figure 80. 

The realization of a l6-point correlator using two CZT ICs is shown in 
Figure 82. This correlator convolves the real input signal with an impulse 
response h(t). True correlation is obtained, since the pulse Tj can blank the 
CCD input when required. 

Another application (not shown) is the realization of a 16-point pro- 
grammable transversal filter. This realization would require eight CZT ICs 
arranged as four correlators of the type shown in Figure 82. T^ and will 
provide synchronization between the chips. These correlators would have the 
inputs tied together and their outputs summed. 

These are only a few of the many applications the general-purpose CZT 1C 
can perform, but they demonstrate its tremendous computational power. 
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SECTION III 


64-point CCD ANALOG/ANALOG CORRELATOR 


A . Design 

The convolver is organized according to the block diagram of Figure 83 . 

Two 64-stage CCDs shift charge in opposite directions. At each stage non- 
destructive outputs are obtained using floating gate amplifiers. Signals are 
buffered, and corresponding outputs from each stage are passed to analog 
multipliers. The outputs of the multipliers are summed to complete the function. 
Feedback linearization from the first buffer output is provided to remove cir- 
cuit non 1 inear i t ies . 


Parallel outputs are obtained using floating gates, as in Figure 84, A 
floating first poly electrode is placed underneath the phase-two clock. The 
floating gate is capacitively coupled to the CCD channel potential. The out- 
put voltage is attenuated approximately 50% due to parasitic capaci tances , 
mostly capacitance to the phase-two clock. The floating gates are periodically 
reset through small series FETs to a reference. 


The analog multiplier is shown in Figure 85 , together with the requ i red 
buffers. Each buffer is biased by a FET current source. The multiplier it- 
self is a two- t rans i s tor circuit, which eliminates the effects of drain resistance 
nonlinearity, as shown below. The current in each of the multiplier FETs as a 
function of the terminal voltages is: 


I 

I 


d2 

dl 


B(V 


gs2 


B(V 


gsl 


't - i'^ds>''ds ^ 


''t - ^''ds’^'ds 


(15) 


where f(V, ) 
ds 

into account 
It is only a 


is the next higher order approximation for the current, which takes 

the transistor threshold voltage variation with substrate voltage. 

function of V. and V . Subtracting the two currents we obtain: 
d s 
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Figure 84 .- Floating Gate Reset Technique. 







06 ) 


d2 


- Id1 = 


■ Vi> 


ds 


since the drain and source voltages are identical for the two devices. Note 
that the drain non 1 inearit ies , as well as the higher order terms, have cancelled. 
The residual nonlinearities that remain arise from variation in FET transcon- 
ductance with gate voltage and from buffer output impedance. 

A feedback linearization circuit has been incorporated into the input stage 
to control the dc voltage levels and eliminate non 1 inear i t ies in the CCD input 
to buffer output transfer function. It is shown in Figure 86. The input stage 
is identical to the other CCD stages and has its own floating gate output and 
buffer. Charge is input to the CCD stage through a FET, where the input signal 
is on the source (the input diode), and the feedback signal is on the gate (the 
feedback gate). As the charge fills up the well in the CCD input stage, the 
floating gate senses it, and the voltage level of the feedback gate drops until 
it is a threshold above the voltage on its source. Thus, the relative input 
and buffer output voltage levels depend only on the characteristics of the input 
FET, and not on all the other adjustments required to make the device operate, 
such as the floating gate reference voltage. 


Note that the charge is metered into the input well in exactly the right 
amount to provide a linear output from the buffer, and the buffer nonlinearity 
caused by body effect and static loading becomes unimportant. 


B . Performance Limitations in Design 
1 . Multiplier Nonlinearity 

There are two main sources of multiplier nonlinearity. The first 
is due to the variation in FET transconductance with gate voltage. It comes 
from carrier velocity saturation effects to a small degree, but mostly from the 
variation in surface mobility as a function of surface concentration. That 
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Figure 86. -CCD Input Stage. 




effect is not generally modelled on SPICE simulation programs accurately, and 
in any case, such simulation was not available to the designer when the device 
was made. The only guidelines available were to attempt to keep the FET gate 
voltages relative to threshold fairly small, such as 3 or 4 V. Distortion from 
this effect was estimated at a fraction of a percent for 1 V signals. 

The second major source of distortion came from the fact that the 
multiplier FETs did not have a voltage source on their drains, but a FET buffer 
with a finite output impedance. This resistance allowed the output voltage to 
be modulated by the loading. This problem was reduced by the feedback lineariza- 
tion circuitry. The buffer of the input stage was loaded with two FETs, both 
with their gates at the reference, or zero voltage level. The effects of this 
"average" load were therefore eliminated. However, the buffer output voltage 
could still be modulated by changes in the loading that occurred whenever a 
signal was presented to the gate of the multiplier FET. This would cause an 
attenuated gate signal to appear at the drain. Whereas the desired output 
current is of the form; 


I = B V V . 
O g a 


V = V o - V , 
g gs2 gsl 


When V, is modulated by V , V , becomes 
d q d 


''d = ''d 


and the output becomes 


I = B V V. (1 - B Z, V ) 
o m g d ' m b g' 


= B (V V . - B Z, V V.) 
m g d m b g d 


(17) 


( 18 ) 


(19) 
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The amount of quadratic distortion is 


D 


A B 


V 


b gst(b) 


( 20 ) 


where A is a factor that accounts for the increase in FET output impedance due 
to the capacitive source impedance at the gate: 


A = (C + C )/C 

source gs source 


( 21 ) 


The amount of distortion can be decreased by making the buffers wide (low 
impedance) and the multiplier FETs high impedance. The dimensions used in the 
design were 7.2 mils wide x 0.3 long for the buffers versus 0.5 W x 2 mils long 
for the resistor FETs. The B-ratio was S6 , resulting in calculated distortion 
of about a percent at 1 V signal levels with 500 pA buffer bias current. 


2 . Offsets 

A severe performance limitation comes from the variations in offset 
voltages of the 64 buffers on each CCD. This adds a random offset term to each 
term of each multiplication, limiting filter coefficient definition, and thus 
filter stopband suppression. It was estimated that the offsets would be about 
30 mV rms , or about 30 dB below the 1 V signals. Assuming the offsets to be 
totally uncorrelated, the processing gain of 8 (square root of N) would drop 
the sidelobe levels to -48 dB. 


However, it was not expected that the offsets would be completely 
uncorrelated , though a design technique was used that helped that to occur. The 
buffer current sources were geometrically identical to the buffers themselves, 
so offsets in the current sources would tend to cancel those of the buffers. 
Since they were within 0.4 mil of each other, reasonable correlation was ex- 
pected between those two offsets. This correlation will reduce the magnitude of 
offsets and decrease the correlation between the offsets of adjacent buffers. 
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3. 


Input Stage Offsets and Non 1 inear it ies 


There is still a residual nonlinearity left in the transfer function 
of the feedback linearized input stage. Since the output equilibrates at a 
threshold above the input, the non 1 inear it ies in threshold voltage versus input 
level would be added to the signal. These were expected to be, however, on the 
order of only a few tenths of a percent. 

Further problems were expected from the input stage in terms of defining 
the absolute dc level of the output. The threshold voltage shift would have to 
be compensated at the input to bring the output to the desired level. This 
level would then have the temperature dependence of a FET threshold voltage. 

In this type of input stage the output voltage approaches the correct level 
asymptotically as 1/T. Assuming that it will never really get there, there 
will be a frequency dependent offset at the buffer outputs. 

C. Measured Performance 

1 . CCD Signal Capacity 

The output voltage swing available at the floating gate buffer out- 
puts was found to be significantly smaller than the design specified. Although 
2 V had been desired, only about 0.8 V was generally measured. Several 
factors were responsible for this severe problem. First, the CCD had been 
designed for 5 \/ implanted wells, but the measured change in threshold voltage 
was only 3.2 V. This would translate to about 2.7 V change in the CCD surface 
potential, or about 85% of the threshold change. 

The interlevel (between first and second poly) oxide thickness was 
reduced between the time the device was designed and the time it was processed. 
This increased the parasitic loading on the floating gate and decreased the 
signal gain. At the calculated gain of 50%, the floating-gate output would be 
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The charge capacity of the phase-one well was smaller than the phase- 
two well that contained the floating gate, partly because the reduced interlevel 
oxide thickness increased the charge capacity of the phase-two well and partly 
because of design error. The charge capacity of the phase-one well was only 
about 80 % of that of the well with the floating gate. That 20% reduction in 
output voltage swing would reduce the available voltage to 1.09 V. 

Finally, with a buffer gain of 0.8, the buffer output swing would be 

0.87 V , which is in agreement with experiment. 

P P 

This limitation can be overcome in a future design, since it was 
imposed mostly by the two-phase structure and by the inappropriate sandwiching 
of the floating gate under another gate, which produced much attenuation of the 
signal . 


In general, the signal levels used were on the order of +/- 1/3 V. 

The measured dynamic range, limited by signal size on the top end and noise on 
the bottom end, was 75 dB, using external BIFET opamp output amplifiers. In 
a future design with 1 V signal levels at each multiplier input, the dynamic 
range would be increased by 20 dB. 

2 . Feedback Linearized Input Stage 

The feedback linearized input stage performed as designed. Figure 87 (a) 
shows a triangle wave at the CCD input with the identical triangle wave offset 
by a threshold at the feedback gate output. The feedback gate waveform is a 
sampled data signal, and the reset portion of each cycle is evident in the 
photo. Note that this signal is about 1.9 V , much larger than the CCD is 
capable of handling. This signal is permitted at the input stage by overfilling 
the well. The barely visible intermediate voltage waveform in the photo shows 
the output voltage after charge transfer when the excess charge is left in the 
input stage, since all of it could not be transferred out. 
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gure 87(a).- CCD Input Stage Feedback-Linearized Outpu 



nsMHiSH 


Figure 87(b).- CCD Serial Output vs Input. 



Figure 87(b) shows a photograph of the CCD input feedback gate wave- 
form, together with the waveform at the serial output buffer of the CCD. Note 
that the waveform is delayed by about 2.5 ms, consistent with the 25 kHz clocking 
frequency of the CCD, and that it saturates in the negative direction, showing 
the limited charge capacity of the CCD. Its output is about 0.8 

The linearity of the transfer function from the CCD input diode to the 

serial output buffer, which should have the same waveforms as the feedback gate, 

was measured using a spectrum analyzer. The particular chip used had a smaller 

signal capacity than observed on chips from other slices. Output voltage signal 

swing was 0.6 V for CCDl (which drove the drain of the multiplier FETs) and 
P P 

0.7 V for CCD2. The spectrum analyzer outputs showed good linearity until the 
PP 

limits of signal swing were reached. At that point the distortion increased 
due to CTE loss in the case of nearly empty wells, or due to a combination of 
CTE loss and CCD flooding in the case of overfull or nearly full wells. 


Figures 88(a) and 88(b) show CCDl output linearity with a 4 V and 

PP 

0.5 Vpp signal swings, respectively. Distortion in the two cases was 0.2% and 

0.28%. Figures 89(a) and 89(b) show CCD2 output linearity with 0,4 and 

0.55 V signal swings. Distortion in the two cases was 0.27% and 1.5%. 

P P 


It is believed that performance of the input stage circuit was good, 
but signal swings were small enough to prohibit truly demanding tests. It is 
also believed that the circuit had enough problems in terms of defining dc 
levels, frequency dependent offsets, and the possibility of instability in some 
applications that its duplication in other circuits is not recommended. Other 
similar circuits will be discussed later. 


3 . Multiplier Linearity 


The multipliers were tested using the circuit of Figure 90. The CCDs 
were disabled, and the buffer FETs were driven through the floating gate reset 
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Figure 88(a).- CCDl 


Serial Output Linearity 0.4 


S igna 1 . 



Figure 88(b). -CCDl Serial Output Linearity 0.5 V 
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Figure 89 (b).- CCD2 Serial 


Output Linearity 0.55 
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Figure 90.- Mu 1 1 i pHer Linearity Test Circuit. 



FETs from the floating gate reference. This circuit was placed inside the 
feedback loop of an opamp, which sampled the serial output of the CCD, to 
accomplish the feedback linearized transfer function. All multipliers were 
tested in parallel. Note that this arrangement allowed the multiplier to be 
tested with larger signals than the CCD would deliver to this circuit. 

The multipliers were tested using a spectrum analyzer. The signal 

on the CCDl buffers (at the multiplier FET drains) was at 2 kHz, while the 

CCD2 buffers were given a 10 kHz sine wave. Figures 91(a) and 91(b) show the 

spectral output with 0.6 V and 2 V sine wave inputs on both terminals. In 

PP PP 

each case, 2 kHz and 10 kHz feedthroughs due to offsets can be seen. The largest 
distortion terms are at 18 and 22 kHz from the buffer output impedance effects. 
Distortion products at 6 and 14 kHz are due to a buffer distortion effect. 


At 0.3 V peak the measured multiplier distortion was about 0.3%, 
while at 1 V peak it measured about 1%. At the 1 V signal levels the multiplier 
reference FET gate voltage was increased to avoid putting the multiplier FETs 
in the saturation region of operation with the larger signals. 


4. Offsets 


The performance of the unit in terms of offset amplitudes was sur- 
prisingly good. Offsets of 30 mV across the entire chip were measured for 

P P 

CCDl. Initially, offsers were measured directly in the time domain. 300 mV 
impulses were sent through both CCDs, giving a correlation peak equal to full 
scale for a single multiplier. This is shown in Figure 92(a), where the correla- 
tion peak has an amplitude of 1 V out of the current differencing amplifiers. 

This waveform also shows the offsets both impulses measure as they shift along 
the chip. 
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Figure 91 (s)- - Mul t ipl ier Output 
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Figure 91 (b).- Mul tipi ier Output = ±1 





Figure 92 (a). - Convol ut ion of Two 300 mV Impulses 
(CCD2 Offsets). 



100 mV/200 M-s 


Figure 92 (b). - Convol ut ion of Impulse with CCD Offsets 
(CCDl Offsets). 
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With impulses sent through only one of the CCDs the measured result 

would be an indication of the offsets of the opposite CCD tap. Figure 92(b) 

shows on a 2X expanded scale the results of those measurements. The top trace 

shows the offsets of CCD2. The offsets are within 60 mV , except for the two 

PP 

endpoints, which are offset outputs from the dummy multipliers at each end of 
the convolver. The bottom trace shows the offsets of CCDl. They are within 


The offset performance of CCDl is better than that of CCD2 for two 
reasons. First, the CCDl output buffers must drive the drains of the multiplier 
FETs and therefore have a rather large geometry. The output buffers of CCD2 
must only drive FET multiplier gates and therefore were physically smaller by 
a factor of 4. Since offset variations would be expected to average out over 
area, the offsets of these smaller buffers should be twice as large as those of 
CCDl . 

The second reason that the offsets of the CCD2 terms were larger is 
that included in those offset terms is any threshold difference between the 
two multiplier FETs. While midway in area between the two buffer types, the 
correspond ing multiplier FET pairs were not as physically close together, which 
would cause threshold variations to have a smaller degree of correlation between 
the two. 


Therefore it is logical that CCDl should have the better offset 
performance . 

Knowing the absolute limits on the offset voltages is not enough to 
describe them adequately, however. Even with a knowledge of the rms value of the 
offsets, it is not known how they would affect the performance of a filter. 

For example, until it is known what the correlation is between the different 
offsets. If it can be assumed that the offsets are uncorrelated, then their 
spectrum is white, and there is a processing gain of the square root of N in 
summing them together. 
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For this reason it is advantageous to make frequency domain measure- 
ments of the offset performance. Such measurements will also give results that 
have a much higher s ignal -to-noise ratio, since the offsets of all the multi- 
pliers are being combined. 

The offsets of CCDl were measured by clocking zeroes through it, 
while running CCD2 from the spectrum analyzer tracking generator. Full scale 
can also be defined by running a large amplitude dc level through CCDl, thus 
implementing a square wave impulse response, or a sin(x)/x frequency response. 
Figure 93(a) shows the resulting outputs with full scale defined as a 0.3 V 
square box response in CCDl. The CCD clocking frequency is 25 kHz. The top 
trace shows the sin(x)/x frequency response, while the bottom trace shows the 
response of the offsets alone. Except below 1 kHz, the response is 50 dB 
below full scale. The increased amplitude at low frequencies is due to very 
slow variations in offset across the length of the chip. Figure 93(b) shows 

the same response out to 10 kHz, and Figure 9^(a) shows the same result out to 

25 kHz. The offset output is at a level of -50 dB out to the Nyquist frequency. 

Note that in Figure 94(a), there are strong offset components at the Nyquist 

frequency. This is due to the buffers being a mirror-image design: alternate 

buffers are mirror images of each other, so that the two power supply buses can 
be shared. This mirror-image design should therefore be avoided in future imple- 
mentations - 

Figure 94(b) shows the same measurements applied to determining the 
offset performance of CCD2. Offsets are about 10 dB worse than for CCDl. 

Figures 95(a) and 95(b) show CCD2 output for full scale and offset output, 
respectively, out to the sampling frequency. 
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-0 dB 


Figure 93(a). -CCDl Frequency Domain Offset Measurements 
(Low Pass Output to 2 . 5 kHz). 



-0 dB 


Figure 93(b). -CCDl Low Pass and Offset Error Outputs to 10 kHz 
(Offset Error Output). 
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-0 dB 



Figure 9^(a).-CCD1 Low Pass and Offset Error Outputs to 2.5 kHz. 



-0 dB 


Figure 9^(b). -CCD2 Low Pass and Offset Error Outputs to 2 , 5 kHz, 
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-0 dB 


Figure 95(a).- CCD2 Low Pass Output to 25 kHz. 



-0 dB 


Figure 95(b). -CCD2 Offset Error Output to 25 kHz. 
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5. 


Multiplier Scale Factor Matching 


Another measure of performance of the multiplier is the degree of 
matching of multiplier scale factors across the length of the chip. This match- 
ing was found to be within about plus or minus 1% over the entire chip. It 
was measured by sending an impulse through CCDl, with either plus or minus full 
scale in CCD2, The difference between these outputs is the multiplier scale 
factor. It was found that the scale factor changed very slowly, almost linearly 
across the chip, as if it were due to a very slow change in the oxide thickness 
in the multiplier FETs. Since the scale factor variations were spatially very 
slow, the test impulse was made four samples wide to increase the S/N ratio. 
Figure 96(a) shows the convolver output when the test impulse was convolved 
with either a plus or minus dc level. The up and down variations of the two 
outputs are due to offsets, while the difference in slopes is due to scale factor 
variations. Figure 96(b) shows the same type of measurement with the dc output 
levels subtracted out. The two outputs should now be superposed. The difference 
in slopes is seen to be very small. 


6 , Eval uat ion 

There were several aspects of performance of this device that were 
interesting and important because they had not been tried before. The CCD 
feedback linearized input stage performed as expected and represents an important 
portion of the system. It also has the potential for use in a modified form in 
other CCD structures. 


The CCD floating gate output performed as designed, although in the 
future it is recommended that alternate layouts be used that do not attenuate 
the CCD signal as much. This attenuation, as well as other design and processing 
problems, limited the signal swings to about 0.8 V , which was the worst 


PP 


limitation of the device. 
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Figure 96 (a). - Mu 1 1 ip 1 ier Scale Factor Accuracy; 

Convolution of CCDl (Impulse of Four 
Samples) with CCD2 (Full Scale Dc). 



Figure 96 (b).- Mul tipi ier Scale Factor Accuracy; 

Multiplier Scale Factor Error. 
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The analog multiplier performed as expected, having about 1% distortion 
level with 1 V signals. The buffer and multiplier offset variations were 
significantly better than expected, and gave perhaps the most important find from 
this circuit: that it is possible to achieve low offset variations in this 

type of correlator, with sidelobe levels typically 60 dB below 1 V signals in 
a 64-stage device. 

Some sample outputs of the convolver in operation are shown in 
Figures 97(a) and 97(b). Figure 97(a) shows the convolution of two triangle 
waves, whose periods are very long compared to the 64 stages of convolution. 

The parabolic outputs are essentially the multiplication of the two triangle 
square waves. Again, the period of the signals is long compared to the convolu- 
tion length, with the result that most of the time the output is multiplying a 
dc level by itself, yielding a positive result. At the square wave transitions 
the transitions approach the center of the convolver from opposite sides; and 
the result becomes the multiplication of two out-of-phase square waves, yielding 
a negative result. 

D . Directions for Future Designs 
1 . System Organization 

It is not felt that there are any advantages at this time to the 
organization wherein both signals are sent through CCDs. More preferable would 
be an organization where the impulse response is loaded onto an array of MOS 
capacitors, which would then drive the multiplier FET drain buffers. This has 
the advantage that a signal can be fed back after the buffer to the sample-and- 
hold circuitry, so that the offset term in the impulse response can be eliminated. 
It would also permit longer storage times of the impulse response. 
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Top and Middle: 
CCDl and CCD2 , 
0.5 V/Div 


Bottom: 

Correlator Out, 
2 V/Div 

5 ms/Div 


Figure 97 (a).- 64-Point Convolution of Two Low-Frequency Triangle Waves, 



Top and Middle; 
CCDl and CCD2 , 
Out, 0.5 V/Div 


Bottom: 

Correlator Out, 
1 V/Div 


Figure 97(b).- 64-Point Convolution of Two Low-Frequency Square Waves. 
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2 . CCD Design 

The design of the CCD will undoubtedly be quite different in future 
devices. The movement at TI to four-phase structures will change the ground 
rules for design radically. The floating gate will probably become one of the 
four phases, guarded on one side from moving clocks by a shield gate. The use 
of four-phase CCDs will result in larger signal capacity, eliminating the small 
signal problem of this device. In any design it will be better not to cover the 
floating gate with another clock electrode. 

3 . Input Stage Design 

It is believed that the use of a linearized input stage is mandatory. 

The type of circuit used here performed adequately, but it would be desirable 
to use a circuit that had no dc level shift from the input to the output, such 
as a potential equilibration CCD input with feedback. 

Care should be taken in the design of the input stage with respect 
to the loading on the output buffer. In this design a dummy stage was added with 
the gates of both multiplier FETs tied to the reference. This will introduce 
signals into the output summing bus unless the thresholds of the dummy multiplier 
FETs are matched exactly. They will not be matched, so the FETs should have their 
outputs connected to a different output ground. 

4. Buffer Design 

To control offsets, it is recommended that large geometry buffers be 
used. This is, of course, necessary in the case of the buffers that drive the 
drains of the multiplier FETs, but should be done also with the other buffers. 

The use of mirror-image geometries from stage-to-stage in the layout 
of the buffers is discouraged because it generates offset components at the 
Nyquist frequency. The buffers and their current sources should be of matching 
geometries and as physically close as possible, so that their offsets will match, 
which will tend to make them cancel. 
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5. 


Mul tipi iers 


Several multiplier topologies are being investigated by various groups. 
They range from a single-transistor, time multiplexed arrangement to a four 
transistor bridge, which has higher linearity. Characteristics of each of 
these circuits are desirable, such as low offsets and high linearity, and the 
different possibilities should be investigated. 

E. Applications and Summary 
1 . Appi icat ions 

It is believed at this time that one of the stronger developing appli- 
cations for ana log/ana log correlators (AAC) is adaptive filtering. It is possible 
that the needed circuitry for updating tap weights could be placed on the chip, re- 
sulting in a device that needed very little external support and could operate with 
high accuracy and dynamic range. The adaptive filtering algorithm would eliminate 
the effects of multiplier offset and scale factor errors and automatically update 
the tap weights. In this application an AAC would be much simpler to use than 
an analog/binary correlator, making maximum use of the strengths of the AAC 
(monotonicity, dynamic range), while eliminating effects from its main drawback — 
offsets . 


While adaptive filtering may be a strong single application, there 
are so many others that the analog/analog correlator should be considered a 
general-purpose part. 

2 . Summa ry 

A 64-stage analog/analog correlator has been designed, built, and 
evaluated. Performance was found to be surprisingly good in spite of a major 
defect that reduced signal voltages by 10 dB from the design goals. The defect 
was in the design of the CCD floating gate outputs and limited signal swings 
to about 0.7 ^pp' Even with those limitations, dynamic range was measured at 
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75 dB, which means 95 dB will be achieved with the proper signal swings. When 
operated as a sin(x)/x lowpass filter, sidelobes due to multiplier offsets 
were 50 dB down from the in-band output. The -60 dB error sidelobes are there- 
fore expected with a corrected design. Multiplier errors were found to be less 
than 1% at the designed signal levels and less than 0.3% at the measured smaller 
levels. The correlator was only operated at low speeds (100 kHz max), where speed 
was entirely limited by the peripheral electronics. 
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SECTION IV 


CCD BINARY/ANALOG CORRELATOR 


A . Introduction 

Classification of mu 1 1 i spectra 1 image data is a complex data processing 
function for satellites that is presently performed on the ground. The goal is 
to perform this pattern classification in the satellite to reduce the data 
transmission rate required for the down link. 

CCD b i na ry/ana 1 og correlators provide a hardware implementation of the 
matrix multiplication function which is the kernel of the pattern classification 
algorithm. At present, the CCD binary/analog correlator appears to offer the 
optimum implementation of the matrix multiplication function. The advantages 
of the binary/analog correlator are given below. 

• High accuracy. Since the weighting coefficients are represented in 
binary form, accuracy of up to 10 bits can be implemented. 

• Arbitrarily long reference retention time. Since the reference is 
stored digitally in static shift registers, it does not need to be refreshed 
as the analog/analog correlator does. 

• The binary/analog correlator is amenable to general-purpose 1C design 
such that a single design can be used in a variety of system applications. 

The problem is to classify vectors X 


( 22 ) 



into classes cu^. . X is an N-element column vector, and 
different wavelength sensors. The value Xj^ represents 
signal in the k^*^ sensor. 


N is the number of 
the amplitude of the 
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A Bayes optimal pattern classifier is one which classifies patterns X into 
classes ujj , j = I, m in such a way as to maximize the discriminant 

g.(X) = - S p(x|cD ) p(u).)/p(X) . ( 23 ) 

jVi -J 


In the above equation 

• p(x|ojj) is the conditional probability density function which gives 
the probability of occurrence of the pattern X, given that X is in fact from 
the class ; 

• p(uu.) is the a priori probability of oDj ; 

• p(X) is the a priori probability of X. 

Maximizing the above discriminant gj(X) is equivalent to minimizing the expected 
loss L (i) = -g.(X), which is given by 

A I 

(f ) = L p (o). |X) . (24) 

jVi 


In other words, L (i) is the sum of the probabilities of classifying X in the 

A 

incorrect class (j i). Equation (23) is obtained by the use of Bayes' rule 

p(X, u)j) = p(xUj) p(ouj) = p(o)J X) p(X) , ( 25 ) 


f rom wh i ch 


p(cu.|x) = p(x|uj ) p(cuJ/p(X) . (26) 

J -i J 

Since p(X) is not a function of i, maximizing g.(X) is equivalent to 
maximiz Ing 

g^(X) = - S p(x|gu,) p(to ) , ( 27 ) 
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which, in turn, can be written as 


g, (X) = p(xlti).) p(cu.) - p(X) . 


Equation (29) is maximum if, and only if, 
g’! = p(x|cuj ) p(u)j ) 

is maximum. Thus, the decision rule is: X e cUj if, and only if, 

p(x|(iJ.) p(u).) ^ p(x|o)j) p(cuj j 7^ i . 

This is commonly referred to as the maximum likelihood decision rule. 


A quadratic classifier utilizes a Gaussian probability distribution 
function. In one dimension this takes the form 


. 2 r (x - Pj) 

p(xlcuj) = (2m ) ^ exp - 2 

' ^ 2a, ^ 


Following Equation (29), it is desirable to maximize the discriminant 


P(«>;) (x - p )' 

9^' (^) = “Z= exp - 2 — 

a, ^ 2a, 

I I 


Maximizing g'/ (x) is equivalent to maximizing 

g"'(x) = log g'/ (x) = log p((u.) - ^ log 2n 


- 1 og Q . - 


(x - Pj) 


I 



Since the constant term in -g log 2tt appears in all of the gf(x), it can be 
dropped, yielding the discriminant 


(x - 

g'J"(x) = log p(cu.) - log a. 2 * (3^) 

2a. 

I 


Returning now to N~d imens iona 1 space and assuming that all classes are 
equally likely on an a priori basis, the discriminant reduces to 


where 


and 


g, (X) = -i logls.l -i(x - U,)"^ S."' (X - u.) 


1 1 

12 

a. 

1 

a . 
1 

21 

22 

a. 

1 

a . 
1 


a 

a 


IN 

2N 


a . 


N1 


N2 


NN 




(35) 


(36) 


(37) 


The simplified discussion given above indicates the importance of matrix 
multiplication to the pattern classification algorithm. In particular, compu- 
tation of the form x"*^CX forms the kernel of the pattern classification algorithm. 
This computation can be implemented with binary/analog correlators. 
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Section IV. B discusses the binary/analog correlator technology and presents 
results obtained with a 32-stage x 4 bit test device. 

Section IV.C discusses the application of the binary/analog correlator to 
the computation of X^CX. 

Section IV.D discusses the preliminary design of an IC to implement a 
classifier having l6 sensors with 8-bit accuracy. 

Section IV.E discusses the design and implementation of a 32-stage by 
8-bit binary/analog correlator unit suitable to implement a classifier having 
M = 8 bits and N = 4 sensors using two of the existing 32-stage by 4-bit devices. 


B . B inary/Ana 1 oq Programmable CCD Correlator Technology 

The concept of the binary/analog correlator is illustrated in Figure 98. 

The weighting coefficients (h^h. . . . h^ ^ of Figure 1 ) are made electronically 
programmable by decomposing each one into a binary representation. h^ is repre- 
sented with M-bit accuracy by 


h 

n 


M-1 

S 

k=0 



n 



(38) 


For M-bit accuracy, M parallel CCD binary/analog correlators are put on a chip 
as shown in Figure 98. The most significant bit h^ of each coefficient is 
loaded into the static shift register (coefficient store) shown as elongated 
rectangles in the filter at the top of the figure. The second most signifi- 
cant bit h^ of each coefficient is loaded into the second coefficient store, 

M-1 

and the least significant bit h is loaded into the coefficient store shown 

at the bottom of the figure. The analog input signal to be filtered is applied 

without attenuation to the top filter (most significant bit). The input is 

attenuated by a factor of two at the input of the second filter (second most 

M-1 

significant bit), and it is attenuated by a factor of 2 at the input to the 
bottom filter (least significant bit). This attenuation is performed using 
capacitive ratio techniques similar to those employed in MOS MDAC technology. 
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Figure 98.“Block Diagram of the Binary/Analog Correlator, 
















The coefficients stored in the static shift registers then control the weight- 
ing of the transversal filters, and when the outputs of each filter are summed . 
together as shown, the result is 


../V '^out^^^ |0 ^-n „-] 1 -n 

in n =0 n =0 


N- 


+ S h”"' z‘" (39) 


n=0 


N-1 ^M-1 


= s r s (h*^ 2"'^y 

^ L,, n V n /. 


n =0 ‘-k^O 


-n 


m 


N-l 

= L h z 

A ^ 

n =0 


-n 


( 41 ) 


In the programmable CCD correlator the relative timing of charge transfer 
is selectively programmed. The structure is sketched in Figure 99- It con- 
sists of a conventional serial ^ CCD shift register in which phase-and-a-ha 1 f 
clocking is employed. All of the electrodes are connected to a common output 
bus line, which is held by the output circuitry to an intermediate potential, i.e,, 
between the OFF potential (0 V) and the ON potential (15 V) . The 02 electrodes 
are individually clocked by a program stored in a digital shift register that 
is physically parallel to the CCD channel. The 0 ^ electrodes are used as barriers 
at the time the charge is sensed and are clocked with a common bus line. The 0 ^^ 
electrodes are biased at a constant potential slightly above 0 V by a common 
bus line. The operation is detailed below. 


Assume that at some time in the clock cycle t^ (see Figure 99) all the 
signal charge packets in the CCD reside under 02 electrodes that are all in 
the ON state. At a later time t 2 some of the 0 ^ electrodes are selectively 
turned OFF (those corresponding to program bits that are equal to 1 ' s) . 
Immediately after t 2 , the signal charge packets will transfer from the OFF 
02 electrodes through the 0 ^ electrodes to the corresponding successive 0 ^ 
electrodes which are maintained at an intermediate potential between the ON 
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Figure 99.- Schemat ic of the Operation of the Programmable CCD Correlator. 



r 


0^ level and the OFF 0^ level. The transferred charges will be sensed and summed 
by the output integrator. Later in the cycle at t^ , the remaining 0^ clocks are 
returned to ground potential and the charge packets that did not transfer at t^ 
will transfer to 0^ electrodes. At time' t^ the 0^ barrier electrodes and 02 
programmed electrodes turn ON, transferring the charges from the 0^ sense 
electrodes. At time t^ the 0^ electrodes turn OFF, leaving the charges under the 
02 electrodes and the cycle is repeated. The 0^ electrodes serve as a buffer to 
prevent the transition of the electrodes at time t 2 from capacitively coupling 
to the integrator inputs (0/^)- Charge is never stored under the 0^ electrodes. 

In the time interval between t 2 and t^ the output signal is available and is 
given by 


^out (n) 


N 

E h X q (n - i ) . 
i=l ' ^ 


( 42 ) 


As described thus far, the device required unipolar signals. To allow 

signal of either sign and to allow the CCD to operate with a fat zero charge 

for better charge transfer efficiency, it is desirable to operate two structures 

like the one just described in parallel in a differential mode. A fat zero 

charge Q_j-^ is inserted into each size of the differential pair, and a signal 

charge Q . is added to one side and substracted from the other side. The 
^ s ig 

output charge will be given by 


^ou t (n ) 


N 

= E 
i = l 


[hi X + 




[h, X (a 


fz 


- 




^ou t (n ) 


N 

E 2h. 

I 


s ig 


( 44 ) 
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For some applications a single-bit coefficient is adequate, but for appli- 
cations requiring coefficients with multiple-bit accuracy, several differential 
channels may be binary weighted and combined as illustrated in Figure 100. This 
figure illustrates a configuration for correlating N analog samples with N 
digital words, each of which is represented by M bits. In this scheme the 
weighting is applied to the analog signals at the input of the CCD shift 
registers. The summing of the M parallel channels is accomplished simply by a 
parallel connection of the clock buses and of the 0^ clock buses. Thus, 
only one output integrator is needed. 

The 02 electrode clocking will be accomplished by the same shift register 
that stores the digital code. The circuit diagram for one stage of the shift 
register is shown in Figure 101. This shift register operates with three non- 
overlapping clocks as shown. Another signal R is used to control whether the 
shift register is connected serially for loading or whether each bit Is con- 
tinuously recirculated within the same stage in a storage mode. 

The device illustrated in Figure 102 has been built to demonstrate the 
b ina ry/ana 1 og concept discussed above. The device has 32 programmable taps, each 
having 4-bit precision. The experimental accuracy of the taps, illustrated in 
Figure 103, is ±1% consistent with LSB of 8-bit taps. The experimental impulse 
response shown in Figure 103 has a digital code -1 -7/8 -6/8 ... +7/8 in order 
to estimate weighting coefficient accuracy. 

Figure 104 illustrates operation of the test device as a bandpass filter. 
Figure 104(a) is a photo of the impulse response and Figure 104(b) is of the 
frequency response. The center frequency is J the clock frequency (f^ = 5 kHz), 
and the impulse response is weighted with a Hamming window quantized to 4-bit 
accuracy. The sidelobes in Figure 104(b) are within 1 dB of the sidelobes 
predicted by a computer simulation of this filter with 4-bit tap weight 
quantization. Figure 105 shows the frequency response of the bina ry/ana log 
correlator as a bandpass filter with rectangular weighting operating at a 500 kHz 
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One Differential Channel 
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Figure 100,- Schematic of a Programmable CCD Correlator with 
M-bit Accuracy on the Weighting Coefficients; 
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Figure 1 02.- Photograph of a 32-Stage Binary/Analog Correlator with 
Four-Bit Weighting Coefficients. This IC has eight CCD 
filters (two for each bit) and static shift registers 
to load and store 32 four-bit words. Two of these ICs 
can be operated in parallel to achieve eight-bit weight 
ing coefficients. Circuitry is included on-chip to 
facilitate microprocessor control. 
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Figure 1 03*- Experimental Impulse Response of 32-Stage 
Binary/Analog Convolver- 
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Figure 1 04.- Response of the Test Device Operated as a Bandpass Filter. 




Figure 1 05.- Frequency Response of the Binary/Analog Correlator 
Programmed to Give a Bandpass Filter with Passband 
at 125 kHz (Clock Frequency is 500 kHz). A rec- 
tangular window function is used. The four bits 
used to code the first five weighting coefficients 
are shown. Note 2 ' s-complement arithmetic is used 
to code the negative values (h ) with b the sign 
bit. ^ 
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clock frequency. Again, the sidelobes are within 1 dB of the computer-simulated 
response. The maximum clock frequency here is limited by the slew rate of the 
external integrating amplifier. 

Figure 1 06 illustrates the use of the correlator for convolving p-n 
sequence codes. The 13“point Barker code impulse response is used. The 3 code 
words used for the p-n weighting are +7/8, 0, and -7/8. Shown are the correlator 
input, the filter impulse response, and the correlator output. The output is 
either 0 or +1 until the correlation peak occurs, at which time the output is 
+ 13. 


C . Application of Binary/Analog Correlator to Pattern Classification 

In a pattern classification system the following operation is needed: 


X^CX 


(45) 


where X is the column matrix of sensor outputs and C = A A is a feature matrix, 

T T T 

Let Y - AX and Y = X A . To evaluate Equation (45) it is sufficient to compute 


where 


T ^ 2 

y'y = (y ) 

i=l 


y ; = E a . . 


X . 


(46) 


(^7) 


i = l 'J J 


The computation of y. can be accomplished using a binary/analog correlator 
described in Section IV. B. 


Figure 107 illustrates the use of the binary/analog correlator structure in 
a system for the computation of Y^Y. The system requires N-stages by M-bit 
binary/analog correlators. N correlators are required for parallel computation 
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Figure 106. 


Matched Filter Operation of the Test 
Device Using a 13-Bit Barker Code. 
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of y.'s. The binary code stored in the i correlator is the i column of 
the A matrix. The output of each correlator is valid after N clock cycles, at 
which time X has been completely loaded into each correlator. The correlator 
outputs are squared and then summed together to obtain Y^Y , 

D . Preliminary Design of Matrix Multiplier IC 

In the desired system, the number of different colored sensors is N = 16, 
and M = 8 bits of accuracy are thought to be required. 

In integrating the system of Figure 107 for 16 sensors, size limitations 
make it desirable to reduce the number of correlators on one chip to four and 
multiplex four chips together. The number of squaring circuits can be reduced 
to one by sampling the correlator outputs after N clock cycles and multiplexing 
these samples off chip. A block diagram of such a chip is shown in Figure 108. 

A preliminary layout of the IC is shown in Figure 109. This IC contains filters 
for four columns of the 16 x l6 feature matrix. The feature coefficients are 
loaded through eight reference inputs, one input for each bit of the reference. 
The reference inputs are TTL levels. A five bit TTL level address (four bits 
for row plus address enable) signal is needed to select the matrix column. An 
input load pulse is needed to transfer the input data to the addressed matrix 
column. The input signal attenuation needed for each column is realized using 
capacitive ratio techniques. The differential current integrator (DCI) consists 
of an MOS amplifier configured as an integrator. The samp 1 e-and-hol d circuitry 
and multiplexer are realized using standard MOS techniques. All clocks are 
derived on-chip from the two-phase master clock. The specification on clock 
rate is 500 kHz (1 MHz master clock). However, the design goal will be 2 MHz 
(4 MHz master clock). Provisions are included for asynchronous operation wherein 
the data can be loaded into the correlators at high speed (greater than 2 MHz), 
and then the clock can be stopped to provide sufficient time 2 \is) to read 
out the matrix product with high accuracy. Synchronizing signals for both inputs 
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l6-*Stage, 8-Bit 
Binary/Analog Correlator 
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Figure 108.- Block Diagram of Binary/Analog Correlator IC. 
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Pin Count 

Feature Coefficients 8 
Analog Input Signal 1 
Master Clock 2 

Address 5 

Output 1 

Sync Pulse 1 

Clear Mode 1 

Voltages: +12 V 1 

-5 V 1 

GND 1 


Total 22 


Figure 109.~ Prel iminary Design of Four-Column by 16-Stage by Eight- 
Bit CCD Matrix Multiplier IC. 


and outputs are provided to allow interfacing to external equipment. Also, 
there is a provision for clearing the stored feature coefficients to zero using 
a digital input control signal. TTL/MOS buffers are provided for all digital 
inputs except the master clock. Sharp transition edges are required, and the 
TTL/MOS buffers required to provide these sharp edges at high speed would 
consume too much power to include on-chip. Dynamic range (maximum output signal/ 
rms output noise) is expected to be 60 dB after the squaring operation. 

Figure 110 illustrates the system configuration. Since this system uses 
parallel processing, an output is available from the accumulator every 16 clock 
cycles. The multiplexer is clocking the l6 correlator samples off chip to the 
squarer and accumulator while the correlators are being loaded with new sensor 
data . 


The matrix multiplier system of Figure 1 08 must be duplicated for each 
feature. The pattern classification decision is then made based on the feature 
matrix multiplier that gives the largest result. 

E . Four-Sensor Demonstration Unit 

Using the 32-stage by 4-bit b ina ry/ana 1 og correlator test devices described 
in Section IV.B, a 32-stage by 8-bit b ina ry/ana 1 og correlator was constructed. 
This unit is suitable to implement a classifier having four sensors and 8-bit 
resolution. The unit is illustrated in Figure 111. To avoid using eight 
chips in this system, parallel processing is not used. All l8 coefficients of 
the 4x4 feature matrix are stored sequentially in the correlator. Two chips 
were required for eight-bit resolution. The system is operated by clocking 
sequences of 32 pulses through the CCDs consisting of the four-term sensor vector 
and 28 "zeros." The accumulator samples every fourth clock cycle. The correla- 
tor output sequence is V f^^ture 1 and V feature 2. A 

squarer and accumulator are required to sum y.^ for i = 1 to 4 for each feature. 
The output will be available once every l6 clock cycles. 
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Figure 111. -Block Diagram of Pattern Classifier Using 
32 X 2-Bit Binary/Analog Correlator. 



Figure 112 illustrates operation of the correlator with eight-bit accuracy. 
Under microprocessor control the digital code is ramped from all +127/128' s 
to all -I's in twos -comp 1 ement form. The output of the correlator (top trace) 
was subtracted from a 10-bit D/A converter with the same input code. The differ- 
ence (bottom trace) was ±20 mV from zero error. Since the binary/analog output 
has a 10 V swing, this represents an accuracy of LSB at eight-bits resolution. 
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Binary/Analog Output 
(10 V Swing) 


Difference Between 
0/4 Converter and 
Binary/Analog Out- 
puts (±20 mV Error) 


Figure 11 2.- 1 1 1 ust rat ion of 8-Bit Accuracy of 
The Binary/Analog Correlator. 
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SECTION V 


CONCLUSIONS 


The objective of this program was to develop different CCD techniques for 
performing programmable correlation for preprocess ing environmental sensor data 
preparatory to its transmission to the ground. Two separate ICs were developed 
under this program and a third IC was evaluated. 

The first IC was a CCD chirp z-transform IC capable of performing a 32-point 
DFT at frequencies to 1 MHz. This IC, suitable for frequency domain correlation, 
contained MOS MDACs, CCD transversal filters, ROM, MOS amplifiers, and all needed 
control circuitry. 

All on-chip circuitry operated as designed with the exception of the limited 
dynamic range caused by a fixed pattern noise due to interactions between the 
digital and analog circuits. Also, the tap weights of the CCD filters were 
erroneously shifted one bit, affecting spectral power density analysis. Both 
problems can easily be corrected with a minor redesign. 

The second IC developed under this program was a 64- stage CCD analog/analog 
correlator for performing time domain correlation. Dynamic range was measured 
at 75 dB and potential for 95 dB dynamic range exists. Multiplier errors were 
found to be less than 1% at designed signal levels and less than 0.3% at the 
measured smaller levels. The main drawback of the ana log/ana log correlator 
is multiplier offsets, but applications exist that eliminate these effects and 
take advantage of the analog/analog correlator's monotonicity and dynamic range. 

During the course of this program, a prototype IC for performing time domain 
correlation was evaluated. This IC offers digital programmab i 1 i ty and appears 
to offer the optimum implementation of the matrix multiplication function needed 
in the pattern classification algorithm. The prototype IC was 32 stages long 
by 4 bits wide and operated at 500 kHz clock rate. The output dynamic range was 
limited to 48 dB by input signal feedthrough, but random output noise was measurec 
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to be 70 dB below the maximum output signal. The feedthrough was caused by 
coupling on chip and can be eliminated in future designs to be 70 dB below 
the maximum output signal. 

It is our strong recommendation that a binary/analog correlator unit be 
developed suitable to implement a classifier having 8-bit accuracy and 16 sensors. 
We also recommend that the applications of the CZT IC to image processing prob- 
lems (two dimensional processing of imagery) be explored. 

In summary, the results of Contract No. NASl-14290 were excellent. The 
first fully integrated CZT IC was designed and evaluated. An ana 1 og/ana 1 og 
correlator was designed and evaluated and an alternate approach to time domain 
correlation, the b ina ry/ana 1 og correlator, was evaluated. A great deal was 
learned under the contract that will enable us to implement, in hardware, a 
pattern classifier suitable for operation in a satellite. 


183 


L 


REFERENCES 


1. Nils J. Nilsson, Learning Machines (McGraw-Hill, New York, 1965). 

2. M. H. White, D. R. Lampe, and J. L. Fagan, "CCD MNOS Devices for Pro- 
grammable Analog Signal Processing and Digital Nonvolatile Memory," 1973 
International Electron Devices Meeting, Technical Digest, pp. 130-133. 

3. J. J. Tiemann, W. E. Engeler, R. D. Baertsch, and D. M. Brown, "Intracell 
Cha rge-Transfer Structures for Signal Processing," IEEE Trans. Electron 
Devices ED-2 1 , 300 (197^). 

4. L. R. Rabiner, R. W. Schafer, and C. M. Rader, "The Chirp Z-Transform 
Algorithm," IEEE Trans, on Audio and Electroacoustics AU- 1 7 > 86 (1969). 

5. H. J. Whitehouse, J. M. Speiser, and R. W. Means, "High Speed Serial 
Access Linear Transform Implementations," presented at the All Applications 
Digital Computer Symposium, Orlando, Florida, January 1973, NUC TN 1026. 

6. R. W. Means, D. D. Buss, and H. J. Whitehouse, "Real Time Discrete Fourier 
Transforms Using Charge Transfer Devices," Proc. CCD Appl. Conf., pp 127- 
139, San Diego, September 1973. 

7. P. Bosshart, "An Integrated Analog Correlator Using Charge-Coupled Devices," 
1976 IEEE International Solid-State Circuits Conf,, Digest of Technical 
Papers, pp 198-199, Philadelphia, February 1976. 


184 


2. Cov«mm«nt Accsttion No. 


3, Rocipionf's Cotolop No. 


4. Titio on4 Swktitio 

LOW POWER, COMPACT CHARGE COUPLED DEVICE 
SIGNAL PROCESSING SYSTEM 


5. Roport Dot* 

February 1980 


6. Porfonntnt Orfoniiotion Co4o 


7. Author<t)p. W. Bosshart, D, D. Buss, W. L. Eversole, 
C. R. Hewes, D. J. Mayer 

9. Porforminp Orponiiotion Homo ond Addrott 

Texas Instruments Incorporated 
Central Research Laboratories 
13500 North Central Expressway 
Dallas, Texas 75265 


12. Sponsorint Agoncy Nomo ond Addrost 

National Aeronautics and Space Administration 
Washington, D.C. 205^6 


I. Porforminf Orfonitotion Roport No. 

08-78-4 


10. Work Unit No. 


II. Centroci or Gront No. 

NAS 1-14290 


13. Typo of Roport ond Ported Covered 

Contractor Report 
April 1976 - August 1978 


14. Sponoerinf Agency Cede 


15. Supplementory Notes 

Langley Technical Monitor: Harry F. Benz 

Final Report 


16. Abstreci The goal of Contract No. NASl-14290 was to develop different charged- 
coupled devices (CCDs) for performing programmable correlation for preprocessing 
environmental sensor data preparatory to its transmission to the ground. Two 
separate ICs were developed under this program, and a third IC was evaluated. 

This final report documents the results of this contract. 

The first IC was a CCD chirp z-transform IC capable of performing a 32-point 
DFT at frequencies to 1 MHz. This IC, suitable for frequency domain correlation, 
contained MOS MDACs, CCD transversal filters, ROM, MOS amplifiers, and all needed 
control circuitry- All on-chip circuitry operated as designed with the exception 
of the limited dynamic range caused by a fixed pattern noise due to interactions 
between the digital and analog circuits. 

The second IC developed under this program was a 64-stage CCD ana 1 og/ana log 
correlator for performing time domain correlation. Multiplier errors were found 
to be less than 1% at designed signal levels and less than 0.3% at the measured 
smaller levels. The main drawback of the ana 1 og/ana 1 og correlator is multiplier 
offsets, but applications exist that eliminate these effects and take advantage of 
the analog/analog correlator's monotonicity and dynamic range. 

During the course of this program, a prototype IC for performing time domain 
correlation was evaluated. This IC offers digital programmability and appears to 
offer the optimum implementation of the matrix multiplication function needed in 
the pattern classification algorithm. The prototype IC was 32 stages long by 4 
bits wide and operated at 500 kHz clock rate. 


17. Kay Words 

Charge Coupled Devices 
Chirp z Transform 
Correlation 


18. Oistribwtion Stotomonf 

Unclassified - Unlimited 

Subject Category 33 


19. Socwrity Clots* f. (of fhis rofrort) 

30 , Socwrity Clotsif. (of this fcfd 

UNCLASSIFIED 

UNCLASSIFIED 



For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley , 1980 


















